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ABSTRACT 


Reflector antennas are widely used in communication satellite systems because 
they provide high gain at low cost. Offset-fed single paraboloids and dual 
reflector offset Cassegrain and Gregorian antennas with multiple focal region 
feeds provide a simple, blockage-free means of forming multiple, shaped and 
isolated beams with low sidelobes. Such antennas are applicable to communica- 
tions satellite frequency reuse systems and earth stations requiring access to 
several satellites. While the single offset paraboloid has been the most 
extensively used configuration for the satellite multiple-beam antenna, the 
trend toward large apertures requiring minimum scanned beam degradation over 
the field of view 18 degrees for full earth coverage from geostationary orbit 
may lead to impractical ly long focal length and large feed arrays. Dual 
reflector antennas offer packaging advantages and more degrees of design 
freedom to improve beam scanning and cross-polarization properties. The 
Cassegrain and Gregorian antennas are the most commonly used dual reflector 
antennas. 

A computer program for calculating the secondary pattern and directivity of a 
generalized dual reflector antenna system has been developed and implemented 
at the NASA Lewis Research Center. The theoretical foundation for this 
program is based on the use of physical optics methodology for describing the 
induced currents on the sub-reflector and main reflector. The resulting 
induced currents on the main reflector are integrated to obtained the antenna 
far-zone electric fields. The computer program is verified with other 
physical optics programs and with measured antenna patterns. The comparison 
shows good agreement in far-field sidelobe reproduction and directivity. 


INTRODUCTION 


The accurate prediction of radiation characteristics for a microwave antenna 
is essential in designing antenna systems. Antenna radiation characteristics 
such as beamwidth, gain, aperture efficiency, side-lobe level, and cross 
polarization are used in analyzing and designing advanced antenna systems. 

The physical optics-current integration approach (ref. 1) described in this 
report is one of several method that can be used for predicting antenna 
performance characteristics. The method assumes that the complex currents in 
both reflectors are known. This, currents satisfies Maxwell's equations and 
are use to solve the complex-vector wave equation at any arbitrary observation 
location. The computation of the induced currents on the main and sub- 



reflector are briefly described. A dual reflector configuration (figure 1.) 
is analyzed, and the results compared with other dual reflector computer pro- 
grams. A description of the input parameters (user guide) and a copy of the 
program are included in Appendixes A,B and C. 


PHYSICAL OPTICS-CURRENT INTEGRATION APPROACH 
DESCRIPTION OF PROBLEM 


The geometry of a dual -reflector with a feed at an arbitrary position is shown 
in Figure 2. Three coordinate systems are shown to define the main reflector, 
the sub-reflector, and the feed position (or array of feeds). The position 
and field vectors of these coordinate system can be interrelated by using the 
Eulerian angles (Figure 3) construction (ref. 2). For instance, the fields of 
the feed can be expressed in feed coordinates (xf,yf,zf) and then transformed 
into sub-reflector coordinates (xs,ys,zs) to determine the scattered field 
from the sub-reflector and then transformed again into main reflector coordi- 
nates (xm,ym,zm) to finally obtain the radiated field of the main reflector. 


INCIDENT ELECTRIC FIELD ON SUB-REFLECTOR 


The radiated electric field from the feed antenna has the asymptotic form 
given by equation (1): 


E( 0,4> ) = F( e,«t> ) 


( 1 ) 


where F( 0,<J> ) is the element pattern, k=2 ti / A is the wavenumber, and r 

is the distance from the source (feed) to the sub-reflector point. The vector 
function in equation (1) can be approximated (ref. 3) by equation (2). 

F(0/<j>) = 0 UE( 0 )(a e^ p cos 4> + b sin <}> ) + 

<)> UH( 0 ) (b cos <(> + a e^ p sin 4> ) (2) 

where UE( 0 ) is the E-plane pattern, UH( 0 ) is the H-plane pattern, and 

a,b, and p are polarization parameters. The various feed polarization 
parameters are described in the following table: 
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TABLE I : Polarization Parameters 


a b p 


Linear x 
Linear y 

Right-hand circular 
Left-hand circular 


1 0 0 

0 1 0 

0.707 0.7007 +90 

0.707 0.707 -90 


Typically these elements patterns can be approximated by a cosine to a power 
function, that is, 


UE( 0 )= cos qe ( 0 ) (3a) 

UH( 0 )= cos^t 0 ) (3b) 

If equations (3a) and (3b) are used to represent the element pattern, the 
power radiated (ref 3.) by this source is given by equation (4). 


Prad 


(ge+gfo+1) 
60(2ge+l) ( 2qh+l ) 


(4) 


SURFACE CURRENT APPROXIMATION 


The foundations of physical optics (PO) , rests on the assumption that the 
induced current on the reflector surface is given (for a perfect conductor)by 

J = 2 (n x H inc ) illuminated region 

J = 0 otherwise 

where n is the unit normal to the surface and H is the incident magnetic 
field. This incident field may emanate directly from the source or be 
scattered from the sub-reflector. Although the PO current is an approximation 
for the true current on the reflector surface, it nevertheless gives very 
accurate results for predicting far-field patterns of reflectors. 
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SCATTERED FIELDS FROM SUB-REFLECTOR 


For a given point on the sub-reflector (xs.ys.zs) and the feed located at 
(xf,yf,zf) the incident fields on the sub-reflector are given 


E = 



F(xs,ys,zs) 


(5a) 


where F(xs,ys,zs) is the feed pattern, k the wavenumber and r the distance 
from the feed to sub-reflector point. The magnetic field incident on the sub- 
reflector is given by 


H = (r x E)/Z 0 (5b) 

The scattered fields from the sub-reflector are given by (ref. 4) 


H(xm,ym,zm) 


= jk ff (J x rl) 


e~ JJcrI 
4 nrl 


ds 


(6a) 


E(xm,ym,zm) 


-jkZ 0 // (J-(J • rl)rl) £—1 ds 


(6b) 


Where J is the induced current on the sub-reflector, rl is the distance from 
any point in the sub-reflector to the observation point (xm,ym,zm). rl is a 
unit vector in the direction from any point in the sub-reflector to any 
observation point on the main reflector (xm,ym,zm). 


MAIN REFLECTOR FAR-FIELDS 


The resulting induced currents produced by the sub-reflector scattering the 
main reflector are integrated to obtain the far-zone electric fields. 
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// (J-(J • R)R) 


ds 


(7a) 


E( 0,4> ) = -jkZ 0 e - jkR 


e jkz 

4 nR 


H( 0,4, ) = (R x E)/Z 0 (7b) 

Where J is the induced current in the main reflector, R is a unit vector from 
any point in the main reflector to the far-field observation point, r is the 
distance from the origin of the main reflector coordinate system to any point 
in the main reflector. 

This method of calculating secondary pattern is accurate in cases where the 
antenna diameter is of the order greater than 50 to 100 wavelength. If the 
antenna diameter is of the order less than 50 wavelength, the accuracy is 
reduced, specifically in the sidelobe region. The reflector configuration 
described in figure 1 was analyzed by using various methods and computer 
codes. The calculated E- and H- plane far-field antenna pattern and directiv- 
ities are shown in figures 4a and 4b respectively. The directivity and the 
far-field pattern are in a very good agreement among computer programs. The 
computer program given in appendix C was used to analyze the configuration. 


DIRECTIVITY 


The far zone electric field is usually divided into two orthogonal polariza- 
tions. Following Ludwig's definition 3 (ref. 4) the following unitary 
polarization vectors are introduced 


R = 0 (a e^P cos 4> + b sin 4> ) + (8a) 

4> (-a e Jp sin 4> + b cos 4> ) 

C = 0 (a e^P sin 4> - b cos 4> ) + (8b) 

4> ( a e^ p cos 4> + b sin 4> ) 

if the secondary pattern can be expressed as 


e = (^(0,4>) +V0,4>) ) 


(9) 
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The reference-polarization expression is 


-ref = E • (**)’ 


(10a) 


and the cross-polarization expression is 


E cross _ E ’ ( E *)* (10b) 

The directivity for the reference polarization is defined by 


DR( 0,4> ) = — 

P rad 

similarly the directivity for the cross polarization is defined by 


(11a) 


DC( 0 , <J> ) = 


^ ^ ^ Cl 033 ^ C l OSS * )/Z, 


Pi ad 


(12b) 


CONCLUDING REMARKS 


A computer program using physical optics-current integration method, has been 
developed for calculating the far-field antenna pattern of dual reflector 
antennas illuminated by a feed with arbitrary polarization. The program 
utilizes a 3th order polynomial spline or nth order polynomial interpolation 
algorithms for cases in which the reflectors are numerically specified. The 
results for the far-field sidelobes and directivity are in good agreement with 
those obtained by other well-known techniques. 

The computer program based on physical optics-current integration techniques 
is one of the main system engineering tools used at NASA Lewis Research Center 
for analyzing advanced antenna systems. 


6 



APPENDIX A 


IDEAL REFLECTOR CONFIGURATIONS 


Offset dual -ref lectors are carved-out of portions of surfaces of revolutions 
(paraboloids t ellipsoids,hyperboloids, etc.) resulting from the intersection 
with cylinders or cones. The cylinders have their axes parallel to the axes 
of the parent reflector surfaces and the cones have their tips at one of the 
foci of the reflectors. In this appendix the geometrical characteristic of 
offset conic sections are presented. 

The following are the analytical equations describing parabolic, hyperbolic 
and elliptical surfaces of revolution all are shown in main reflector coordi- 
nate system. 


A: Parabolic reflector : The geometry associated with a parabolic 

reflector is shown in figure A-l 


Parameters : F focal length 


Equation 


x 2 +y 2 

4F 


B: Hyperbolic Sub-reflector : The geometry associated with a hyperbolic 

sub-reflector is shown in figure A-2. 


Parameters 


Zq offset distance 

a vertex distance 

b = fc^a 2 
2c foci distance 


Equation 


z = 




1 + 


x 2 +y 2 

b 2 
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Elliptical sub-reflector : The geometry associated with an elliptical 

sub-reflector is shown in figure A-3. 


Parameters 


Equation 


Zq offset distance 
a vertex distance 


b = \/a 2 -c 2 

2c foci distance 


z = 


z 0 +a, 


N 


_ x 2 +y 2 

b 2 
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APPENDIX B 


PROGRAM INPUT USER GUIDE 


A computer program was designed to calculate the antenna performance charac- 
teristics. The method of analysis is physical optics. The program runs in an 
IBM370 using VM operating system. All the inputs are put into the program 
DRSG FORTRAN and are describe as follows: 


FFREQ 

QQ 

DMX.DMY 

DSX.DSY 

MAXMX.MAXSY 
xmO,ymO,zmO 
xsO,ysO,zsO 
xf ,yf ,zf 
xr.yr.zr 
rtempl (sub) 


frequency GHz 

feed pattern exponent 

x and y length in wavelength main reflector 
rectangle 

x and y length in wavelength sub-reflector 
rectangle 

number of points in the x and y direction 
lower left corner of main reflector rectangle 
lower left corner of sub-reflector rectangle 
feed location in wavelength 
feed boresight point on sub-reflector 
parameter a in wavelength 


rtemp2 

rtemp6 
fradsq 
FCENX, FCENY 
radsq 

CNTRX ( CNTRY 
rtempl (main) 


parameter J b= v / a 2 -c 2 

offset distance in wavelength 
radius square of cylinder sub-reflector 
center of sub-reflector cylinder 
radius of cylinder of main reflector 
center of cylinder of main reflector 
1/4F, F is focal length in wavelength 
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APPtNDIX C 
COMPUTER PROGRAM 


PROGRAM DRSG 

******************************************************************** 

* NASA LEWIS RESEARCH CENTER * 

* AUTHOR : R. ACOSTA AND A. LAGIN IBM 370 VM VERSION * 

* DATE : 7/14/91 * 

* PURPOSE : TO SIMULATE GENERALIZED SURFACES OF REVOLUTION * 

* FOR THE GENERALIZED DUAL REFLECTOR ANALYSIS PROGRAM * 

******************************************************************** 

** COMPILE CONSTANTS 

******************************************************************** 
REAL C 

PARAMETER (C = 2.99792458E+08 ) 

REAL PI 

PARAMETER (PI = 3.141592653589793238) 

REAL ETA 

PARAMETER (ETA = 4.0*PI*C* (1 . E-7) ) 

******************************************************************* 

****** INPUT : FREQUENCY (HZ)************************************** 
******************************************************************* 

REAL FREQQ 

PARAMETER (FREQQ = 19.45*(1 . E+9) ) 

******************************************************************* 

REAL LAMBDA 

PARAMETER (LAMBDA = (C/FREQQ) *1 .000 ) 

REAL LAMBSQ 

PARAMETER (LAMBSQ = LAMBDA*LAMBDA ) 

****************************************************************** 

****** INPUT : QQ : FEED PATTERN EXPONENT************************** 
****************************************************************** 


REAL QQ 

PARAMETER (QQ =62.00 ) 

***************************************************************** 

INTEGER MXYZ 

PARAMETER (MXYZ = 4 ) 

INTEGER DXYZ 

PARAMETER (DXYZ = 3 ) 

INTEGER LNRY 

PARAMETER (LNRY = 1 ) 

INTEGER MDTX 

PARAMETER (MDTX = 10 ) 


****************************************************************** 
** INTRINSIC FUNCTIONS 

****************************************************************** 

INTRINSIC SQRT 
INTRINSIC INT 
INTRINSIC NINT 

****************************************************************** 
** EXTERNAL SUBROUTINES 

****************************************************************** 
EXTERNAL SUBMAI 


DRS00010 
DRS00020 
DRS00030 
DRS00040 
DRS00050 
DRS00060 
DRS00070 
DRS00080 
DRS00090 
DRS00100 
DRS001 10 
DRS00120 
DRS00130 
DRS00140 
DRS00150 
DRS00160 
DRS00170 
DRS00180 
DRS00190 
DRS00200 
DRS00210 
DRS00220 
DRS00230 
DRS00240 
DRS00250 
DRS00260 
DRS00270 
DRS00280 
DRS00290 
DRS00300 
DRS00310 
DRS00320 
DRS00330 
DRS00340 
DRS00350 
DRS00360 
DRS00370 
DRS00380 
DRS00390 
DRS00400 
DRS00410 
DRS00420 
DRS00430 
DRS00440 
DRS00450 
DRS00460 
DRS00470 
DRS00480 
DRS00490 
DRS00500 


******************************************************************** ****qrs005 10 


REAL 

DATARY (MDTX) 

DRS00520 

REAL 

DMX.DMY 

DRS00530 

REAL 

INCMX, INCMY 

DRS00540 

REAL 

DSX.DSY 

DRS00550 

REAL 

INCSX, INCSY 

DRS00560 

REAL 

FREQ 

DRS00570 

REAL 

Q 

DRS00580 


********************************** **************************************QR50O5gO 


** EQUIVALENCE 


DRS00600 
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************************************************************************DR500610 


EQUIVALENCE (DATARY ( 1) 
EQUIVALENCE (DATARY ( 2) 
EQUIVALENCE (DATARY ( 3) 
EQUIVALENCE (DATARY ( 4) 
EQUIVALENCE (DATARY ( 5) 
EQUIVALENCE (DATARY ( 6) 
EQUIVALENCE (DATARY ( 7) 
EQUIVALENCE (DATARY ( 8) 
EQUIVALENCE (DATARY ( 9) 
EQUIVALENCE (DATARY (10) 
INTEGER MAXMX, MAXMY 
INTEGER MAXSX, MAXSY 


,DMX ) 
,DMY ) 
, INCMX) 
, INCMY) 
,DSX ) 
,DSY ) 
, INCSX) 
, INCSY) 
i FREQ ) 

.Q ) 


******************************************************************** 


** INITIAL CALCUALTIONS 

******************************************************************** 


*INPUT : LENGTH OF MAIN REF. RECTANGLE GRID 
******************************************************************** 

DMX = (220.0 * LAMBDA) 

DMY = (220.0 * LAMBDA) 

******************************************************************** 


****** IN p UT . number POINTS IN X AND Y MAIN REF. RECTANGLE GRID***** 
******************************************************************** 


MAXMX =101 
MAXMY = 101 

******************************************************************** 


INCMX = DMX/ (REAL( (MAXMX-1) ) ) 
INCMY = DMY/ (REAL( (MAXMY-1) ) ) 


******************************************************************** 

****** I Np UT : LENGHT OF SUB REF. RECTANGLE GRID ********** 

******************************************************************** 


* LAMBDA) 
* LAMBDA) 


DSX = ( 82.50 

DSY = ( 82.50 

******************************************************************** 

****** INPUT . NUMBER OF POINTS IN X AND Y IN THE SUB REFLECTOR GRID* 
******************************************************************** 


MAXSX = 61 
MAXSY = 61 

******************************************************************** 

INCSX = (DSX/REAL(MAXSX-1) ) 

INCSY = (DSY/REAL(MAXSY-1)) 

******************************************************************** 


FREQ = FREQQ 

Q ■ QQ 

CALL SUBMAI (MAXMY , MAXMX , MAXSY , MAXSX , MDTX , DATARY) 
END 


******************************************************************** 

********* ^SUBROUTINE SUBMA I ***************************************** 
******************************************************************** 


SUBROUTINE SUBMAI ( PMXMY , PMXMX , PMXSY , PMXSX , PDTX , DTAARX) 

INTEGER PMXMY 

INTEGER PMXMX 

INTEGER PMXSY 

INTEGER PMXSX 

INTEGER PDTX 

REAL DTAARX (PDTX) 

REAL C 

PARAMETER (C = 2.99792458E+08 ) 

REAL PI 

PARAMETER (PI = 3.141592653589793238) 


DRS00620 

DRS00630 

DRS00640 

DRS00650 

DRS00660 

DRS00670 

DRS00680 

DRS00690 

DRS00700 

DRS00710 

DRS00720 

DRS00730 

DRS00740 

DRS00750 

DRS00760 

DRS00770 

DRS00780 

DRS00790 

DRS00800 

DRS00810 

DRS00820 

DRS00830 

DRS00840 

DRS00850 

DRS00860 

DRS00870 

DRS00880 

DRS00890 

DRS00900 

DRS00910 

DRS00920 

DRS00930 

DRS00940 

DRS00950 

DRS00960 

DRS00970 

DRS00980 

DRS00990 

DRS01000 

DRS01010 

DRS01020 

DRS01030 

DRS01040 

DRS01050 

DRS01060 

DRS01070 

DRS01080 

DRS01090 

DRS01100 

DRS01110 

DRS01120 

DRS01130 

DRS01140 

DRS01150 

DRS01160 

DRS01170 

DRS01180 

DRS01190 

DRS01200 



REAL 

ETA 



DRS01210 

PARAMETER (ETA 

= 4.0*PI*C*(1 .E-7) 

) 

DRS01220 

******************************************************************* 

DRS01230 

****** input : FREQUENCY******************************************** 

DRS01240 

******************************************************************* 

DRS01250 

REAL 

FREQQ 



DRS01260 

PARAMETER (FREQQ 

= 19.45*(l.E+9) 

) 

DRS01270 

****************************************************************** 

DRS01280 

REAL 

LAMBDA 



DRS01290 

PARAMETER (LAMBDA 

= (C/FREQQ) *1 . OOO 

) 

DRS01300 

REAL 

LAMBSQ 



DRS01310 

PARAMETER (LAMBSQ 

= LAMBDA* LAMBDA 

) 

DRS01320 

****************************************************************** 

DRS01330 

****** input : QQ: EXPONED EXPONENT********************************* 

DRS01340 

******************************************************************* 

DRS01350 

REAL 

QQ 



DRS01360 

PARAMETER 

(QQ 

= 62.000 

) 

DRS01370 

****************************** * **** ***************** * ************** 

DRS01380 

INTEGER 

MXYZ 



DRS01390 

PARAMETER 

(MXYZ 

= 4 

) 

DRS01400 

INTEGER 

DXYZ 



DRS01410 

PARAMETER 

(DXYZ 

= 3 

) 

DRS01420 

INTEGER 

LNRY 



DRS01430 

PARAMETER 

(LNRY 

= 1 

) 

DRS01440 

INTEGER 

VDTX 



DRS01450 

PARAMETER 

(VDTX 

= 4 

) 

DRS01460 

INTEGER 

MDTX 



DRS01470 

PARAMETER 

(MDTX 

= 10 

) 

DRS01480 

INTEGER 

IDXZ 



DRS01490 

INTEGER 

IDXZX 



DRS01500 

INTEGER 

IDXZY 



DRS01510 

INTEGER 

IDXMSK 



DRS01520 

INTEGER 

MXMX, MXMY, TMXMX, TMXMY 


DRS01530 

INTEGER 

MXSX , MXSY , TMXSX , TMXS Y 


DRS01540 

PARAMETER 

(IDXZ 

= 1 

) 

DRS01550 

PARAMETER 

(IDXZX 

= 2 

) 

DRS01560 

PARAMETER 

(IDXZY 

= 3 

) 

DRS01570 

PARAMETER 

(IDXMSK 

= 4 

) 

DRS01580 

******************************************************************* 

DRS01590 

****** INPUT : NUMBER OF 

POINTS (-1) MAIN AND 

SUB REFLECTOR GRIDS*** 

DRS01600 

******************* ************************************************ 

DRS01610 

PARAMETER (MXMX 

= 101 

) 

DRS01620 

PARAMETER (MXMY 

= 101 

) 

DRS01630 

PARAMETER (MXSX 

= 61 

) 

DRS01640 

PARAMETER (MXSY 

= 61 

) 

DRS01650 

******** *********************************************************** 

DRS01660 

PARAMETER 

(TMXMX 

= MXMX + 1 

) 

DRS01670 

PARAMETER 

(TMXMY 

= MXMY + 1 

) 

DRS01680 

PARAMETER 

(TMXSX 

= MXSX + 1 

) 

DRS01690 

PARAMETER 

(TMXSY 

= MXSY + 1 

) 

DRS01700 

************ *********************************************************** 

*DRS01710 

************************************************************************DP^q27 2q 

REAL 

MAIARY(4, TMXMY, TMXMX) 


DRS01730 

REAL 

SUBARY (4, TMXSY, TMXSX) 


DRS01740 

REAL 

XYZARY (DXYZ, MXYZ) 


DRS01750 

REAL 

MRX YZO (DXYZ) , XMO , YMO , ZMO 


DRS01760 

REAL 

SRXYZO(DXYZ) ,XSO, YSO.ZSO 


DRS01770 

REAL 

FEDXYZ (DXYZ) , XF , YF ,ZF 


DRS01780 

REAL 

REFXYZ (DXYZ) , XR , YR ,ZR 


DRS01790 

REAL 

DMX , DMY 


DRS01800 
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REAL INCMX, INCMY DRS01810 

REAL DSX.DSY DRS01820 

REAL INCSX, INCSY DRS01830 

REAL FREQ DRS01840 

REAL Q DRS01850 

************************************************************************DRS01860 

** EQUIVALENCE DRS01870 

************************************************************************DRS01880 

EQUIVALENCE (XYZARY(1 , 1) ,MRXYZ0(1) ) DRS01890 

EQUIVALENCE (XYZARY(1,2) ,SRXYZ0(1)) DRS01900 

EQUIVALENCE (XYZARY(1,3) , FEDXYZ (1) ) DRS01910 

EQUIVALENCE (XYZARY(1 ,4) , REFXYZ (1) ) DRS01920 

EQUIVALENCE (MRXYZO(l) ,XMO ) DRS01930 

EQUIVALENCE (MRXYZO(2) , YMO ) DRS01940 

EQUIVALENCE (MRXYZO(3) , ZMO ) DRS01950 

EQUIVALENCE (SRXYZO(l) ,XSO ) DRS01960 

EQUIVALENCE (SRXYZO(2) , YSO ) DRS01970 

EQUIVALENCE (SRXYZO(3) ,ZSO ) DRS01980 

EQUIVALENCE (FEDXYZ(l) ,XF ) DRS01990 

EQUIVALENCE (FEDXYZ(2) , YF ) DRS02000 

EQUIVALENCE (FEDXYZ (3) , ZF ) DRS02010 

EQUIVALENCE (REFXYZ(l) ,XR ) DRS02020 

EQUIVALENCE (REFXYZ (2) , YR ) DRS02030 

EQUIVALENCE (REFXYZ (3) , ZR ) DRS02040 




a****************************************************************** *****DRS02060 


REAL 

MAIXYZ(3),XM, 

YM, 

ZM 

DRS02070 

REAL 

SUBXYZ(3) , XS , 

YS, 

ZS 

DRS02080 

REAL 

GENXYZ (3) , XI , 

YJ, 

ZIJ 

DRS02090 

REAL 

RADSQ,XSQ 



DRS02100 

REAL 

CNTRX 



DRS021 10 

REAL 

CNTRY 



DRS02120 

REAL 

RTEMPO 



DRS02130 

REAL 

RTEMP1 



DRS02140 

REAL 

RTEMP2 



DRS02150 

REAL 

RTEMP3 



DRS02160 

REAL 

RTEMP4 



DRS02170 

REAL 

RTEMP5 



DRS02180 

REAL 

RTEMP6 



DRS02190 

INTEGER 

I.J.V.W 



DRS02200 


************************************************************************^§02210 
** EQUIVALENCE DRS02220 


EQUIVALENCE (SUBXYZ(l) , XS ) DRS02240 
EQUIVALENCE (SUBXYZ(2) , YS ) DRS02250 
EQUIVALENCE (SUBXYZ(3) ,ZS ) DRS02260 
EQUIVALENCE (MAIXYZ(l) ,XM ) DRS02270 
EQUIVALENCE (MAIXYZ (2) , YM ) DRS02280 
EQUIVALENCE (MAIXYZ(3) ,ZH ) DRS02290 
EQUIVALENCE (GENXYZ(l) ,XI ) DRS02300 
EQUIVALENCE (GENXYZ (2) , YJ ) DRS02310 
EQUIVALENCE (GENXYZ (3) , ZIJ ) DRS02320 


*★*★★***★**★*****★**★*****★*★★★ *****************************************QRS023 30 


** Initialize Arrays DRS02340 
** MAIFIL <== MAIARY () DRS02350 
** SUBFIL <== SUBARY () DRS02360 
** XYZFIL <== XYZARYQ <== MRXYZO() ,SRXYZO() , FEDXYZQ ,REFXYZ() DRS02370 


** DTAFIL <== DTAARXQ <== DMX, DMY, INCMX, INCMY, DSX,DSY, INCSX, INCSY, FREQ, DRS02380 
******************** ****************************************************0RSO239O 


DO 20200 I = 1 , TMXMX, 1 


DRS02400 
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DO 20100 J = l.TMXMY, 1 DRS02410 

DO 20000 W = 1,4,1 DRS02420 

MAIARY (W, J , I) = 0. DRS02430 

20000 CONTINUE DRS02440 

20100 CONTINUE DRS02450 

20200 CONTINUE ' DRS02460 

DO 20500 I = 1 , TMXSX, 1 DRS02470 

DO 20400 J = l.TMXSY, 1 DRS02480 

DO 20300 W = 1,4,1 ; DRS02490 

SUBARY (W, J, I) = 0. DRS02500 

20300 CONTINUE DRS02510 

20400 CONTINUE DRS02520 

20500 CONTINUE DRS02530 


********************************* ***************************************QpsQ2f}4() 
****** jtypyj ; LOWER LEFT CORNER OF MAIN REF GRID************************DRS02550 


************************************************************************ QP5Q2560 
XMO = (-110.00000000 * LAMBDA) DRS02570 

YMO = ( 79.65000000 * LAMBDA) DRS02580 

ZMO = ( 1.000000000 * LAMBDA) DRS02590 

************************************************************************0^Q25oo 

******INPUT : LOWER LEFT CORNER OF SUB REFLECTOR GRID*********************DRS02610 


***************************************************************** 
XSO = ( -41.25000000 * LAMBDA) 

YSO = ( 000.00000000 * LAMBDA) 

ZSO = ( 1.00000000 * LAMBDA) 

*********** ****************************************************** 

****** input : FEED COORDINATES*********************************** 
***************************************************************** 


****** B0RSIGHT BEAM FEED COORDINATES: 
XF = ( 0.000000000 * LAMBDA) 

* YF = ( 0.000000000 * LAMBDA) 

* ZF = ( 103.93560000 * LAMBDA) 

***** CLEVELAND BEAM FEED COORDINATES: 

YF = ( 5.629 * LAMBDA) 

ZF = ( 101.067 * LAMBDA) 

*****MIAMI BEAM FEED COORDINATES: 

* XF = ( -16.4 * LAMBDA) 

* YF = ( 9.62287 * LAMBDA) 

* ZF = ( 99.0325 * LAMBDA) 


*****Lqs ANGELES BEAM FEED COORDINATES: 


★ 

XF = (-7.996 

*LAMBDA) 

* 

YF = (-30.52 

*LAMBDA) 

★ 

ZF = ( 119.49 

*LAMBDA) 

***** 

SEATLE BEAM FEED 

COORDINATES: 

* 

XF = ( 3.4579 

*LAMBDA) 

* 

YF = (-30.2 

*LAMBDA) 

* 

ZF = (119.32 

*LAMBDA) 


************************ **************************************** 


****** INPUT : REFERENCE RAY LOCATION COORDINATES ******* 

**************************************************************** 

XR = ( 0.000000000 *LAMBDA) 

YR = ( 41.25000000 *LAMBDA) 

* YR = ( 30 .O0OOOOOO *LAMBDA) 

Y = YR*39.36 

ZR = ( 11.5*SQRT(l+(XR**2+Y**2)/1058)+97.5) 

ZR=ZR/39. 36 

★**********★★**★*★**★★***★*****★★*******★****★*★**************** 
DMX = DTAARX( 1) 

DMY = DTAARX( 2) 

INCMX = DTAARX( 3) 


DRS02620 

DRS02630 

DRS02640 

DRS02650 

DRS02660 

DRS02670 

DRS02680 

DRS02690 

DRS02700 

DRS02710 


DRS02720 

DRS02730 

DRS02740 

DRS02750 

DRS02760 


DRS02770 

DRS02780 

DRS02790 

DRS02800 

DRS02810 
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INCMY 

= DTAARX( 4) 

DSX 

= DTAARX( 5) 

DSY 

= DTAARX( 6) 

INCSX 

= DTAARX( 7) 

INCSY 

= DTAARX( 8) 

FREQ 

= DTAARX( 9) 

Q 

= DTAARX(IO) 


**************************************************************** 

*******jq G ENEARATE THE INPUT ARRA YS**************************** 
**************************************************************** 

DO 10 1=1,10 
WRITE(15,805)DTAARX(I) 

805 FORMAT (5X, E15.8) 

10 CONTINUE 

DO 11 1=1,4 
DO 12 J-1,3 

WRITE(16,806) XYZARY (J , I) 

806 FORMAT (5X, E15 .8) 

12 CONTINUE 

11 CONTINUE 

************************************************************** 

** Calculate Z, Zx, Zy, Usage for SubReflector. 

****Z = SUBARY (IDXZ, SURFACE Z 

****DX= SUBARY (IDXZX. . . DERIVATIVE WITH RESP. TO X 

****DY= SUBARY (IDXZY DERIVATIVE WITH RESP. TO Y 

************************************************************* 

RTEMPO = ( LAMBDA) 

************************************************************** 

****** jNP UT . A ; PARAMETER OF SURFACE OF REVOLUTION*********** 
*★★************★**★**★*★★********★*★★**********★********★*★★*★****** 

RTEMP1 = (18.97*RTEMP0) 

*******************************************************^****** 

****** INPUT . !/b**2 PARAMETER OF SURFACE OF REVOLUTION******* 
************************************************************** 

RTEMP2 = (1 . / ( .6829) ) 

*************************************************************** 

****** INPUT ; zo OFFSET OF THE CENTER OF THE SURFACE OFF REV.** 
*************************************************************** 

RTEMP6 = 160.85*RTEMP0 

******************************************************************** 

****INPUT: RADIUS OF CYLINDER OF SUB-REFLECTOR********************** 
****************************** ************************************** 

FRADSQ= (19.0*RTEMPO)**2 

******************************************************************** 

****INPUT: CENTER OF COORDINATES OF CYLINDER*********************** 

******************************************************************** 

FCENX=0.0 

FCENY=31.8*RTEMP0 

******************************************************************** 

XS = XSO - INCSX 

DO 20700 I = l.TMXSX.l 

XS = XS + INCSX 
XSQ = XS*XS 
YS = YSO - INCSY 
DO 20600 J = 1 , TMXSY, 1 
YS = YS + INCSY 

FTEMP= (XS-FCENX)**2+(YS-FCENY) **2 
* IF(FTEMP.GT.FRADSQ) GO TO 309 

RTEMP3 = ((XSQ) + (YS*YS)) 

RTEMP4 = (SQRT (1+ RTEMP3*RTEMP2) ) 


DRS02820 

DRS02830 

DRS02840 

DRS02850 

DRS02860 

DRS02870 

DRS02880 

DRS02890 

DRS02900 

DRS02910 

DRS02920 

DRS02930 

DRS02940 

DRS02950 

DRS02960 

DRS02970 

DRS02980 

DRS02990 

DRS03000 

DRS03010 

DRS03020 

DRS03030 

DRS03040 

DRS03050 

DRS03060 

DRS03070 

DRS03080 

DRS03090 

DRS03100 

DRS03110 

DRS03120 

DRS03130 

DRS03140 

DRS03150 

DRS03160 

DRS03170 

DRS03180 

DRS03190 

DRS03200 

DRS03210 

DRS03220 

DRS03230 

DRS03240 

DRS03250 

DRS03260 

DRS03270 

DRS03280 

DRS03290 

DRS03300 

DRS03310 

DRS03320 

DRS03330 

DRS03340 

DRS03350 

DRS03360 

DRS03370 

DRS03380 

DRS03390 

DRS03400 

DRS03410 



RTEMP5 = ( (RTEMP1*RTEMP2) /RTEMP4) DRS03420 

DRS03430 

DRS03450 

| SubRef lector: Hyperbola DRS03460 

DRS03470 

| + DRS03480 

DRS03490 

| / 2 2 DRS03500 

DRS03510 

| Z = 18.97*LAMBDA X / 1 + (X + Y) ) +160.85*LAMBDADRS03520 

DRS03530 
DRS03540 
DRS03550 
DRS03560 
DRS03570 
DRS03580 
DRS03590 
DRS03600 
DRS03610 
DRS03620 
DRS03630 
DRS03640 
DRS03650 
DRS03660 
DRS03670 
DRS03680 
DRS03690 
DRS03700 
DRS03710 
DRS03720 
DRS03730 
DRS03740 
DRS03750 
DRS03760 
DRS03770 
DRS03780 
DRS03790 
DRS03800 
DRS03810 
DRS03820 
DRS03830 
DRS03840 
DRS03850 


DRS03870 

DRS03880 

* DRS03890 

SUBARY (IDXZ ,J,I) = ( (RTEMP6) +RTEMP1*RTEMP4) DRS03900 

SUBARY (IDXZX ,J,I) = (XS*RTEMP5 ) DRS03910 

SUBARY (IDXZY ,J,I) = (YS*RTEMP5 ) DRS03920 

SUBARY (IDXMSK.J, I) = 1. DRS03930 

* DRS03940 

20520 FORMAT (5X , 4(E15 .8, 2X) ) DRS03950 

309 WRITE(17, 20520) SUBARY(IDXZ ,J, I) .SUBARY (IDXZX DRS03960 

1 SUBARY(IDXZY , J, I) .SUBARY (IDXMSK.J, I) DRS03970 

20600 CONTINUE DRS03980 

20700 CONTINUE DRS03990 


DRS04000 

************************************************************************^$ 0 ^^ 


* 

* 

* 

* 

* 

* 

★ 

* 

* 

* 

* 

★ 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

★ 

* 

* 


\/ 


( 0.6829 ) 


+ 

/ 2 2 

X / X + Y 

ZX = 18.97 ( ) / / l+( ) 

150.X LAMBDA \/ 0.6829 

+ 

/ 2 2 

Y / X + Y 

ZY = 18.97 ( ) / / l+( ) 

150.XLAMBDA \/ 0.6829 
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** Calculate Z, Zx, Zy, Usage for MalnReflector. DRS04020 
****Z = MAIARY(IDXZ, .... SURFACE FUNCTION DRS04030 
****Z = MAIARY (IDXZX,... DERIVATIVE WITH RESP. X DRS04040 
****Z = MAIARY (IDXZY,... DERIVATIVE WITH RESP. Y DRS04050 
*********** *************************************************************qrs04060 


RTEMP0= LAMBDA 
XM = XMO - INCMX 

********************************************************************* 

*****INPUT: RADIUS OF CYLINDER MAIN REFLECTOR************************ 
********************************************************************* 

* RADSQ = (107 . 2*RTEMP0) **2 
RADSQ = (108.2*RTEMP0)**2 

********************************************************************* 
****INPUT : CENTER OF COORDINATES OF CYLINDER************************* 

★★a********************************* ********************************* 


CNTRX =0.0 

CNTRY = 189.7*RTEMP0 

********************************************************************* 

*****j^pyj . i/4p PARAMETER FOCAL LENGTH***************************** 
********************************************************************* 

RTEMP1 = (l./(870.00*RTEMP0)) 

RTEMP2 = 2*RTEMP1 
DO 20900 I = 1 , TMXMX, 1 
XM = XM + INCMX 

XSQ = XM*XM 

YM = YMO - INCMY 
DO 20800 J = 1 , TMXMY, 1 
YM=YM+INCMY 

RTEMP4 = ( ( (XM-CNTRX)* (XM-CNTRX) )+((YM-CNTRY)*(YM-CNTRY) ) ) 
IF (RTEMP4 .GT. RADSQ) GO TO 450 


* 

* 


Mai Reflector: Parabola 


* 


* 


2 2 
( x + Y ) 


(870. X LAMBDA) 


* 


* 


X 


Zx = 


(435. X LAMBDA) 


Y 


* 

* 


zy = 


* | (435. X LAMBDA) 

* 


DRS04070 

DRS04080 

DRS04090 

DRS04100 

DRS04110 

DRS04120 

DRS04130 

DRS04140 

DRS04150 

DRS04160 

DRS04170 

DRS04180 

DRS04190 

DRS04200 

DRS04210 

DRS04220 

DRS04230 

DRS04240 

DRS04250 

DRS04260 

DRS04270 

DRS04280 

DRS04290 

DRS04300 

=DRS04310 

DRS04320 

DRS04330 

DRS04340 

DRS04350 

DRS04360 

DRS04370 

DRS04380 

DRS04390 

DRS04400 

DRS04410 

DRS04420 

DRS04430 

DRS04440 

DRS04450 

DRS04460 

DRS04470 

DRS04480 

DRS04490 

DRS04500 

DRS04510 

DRS04520 

DRS04530 

DRS04540 

DRS04550 

DRS04560 

DRS04570 

DRS04580 

DRS04590 

DRS04600 
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MAIARY (IDXZ ,J,I) = ((XM*XM)+(YM*YM))*RTEMP1 

MAIARY(IDXZX ,J,I) = (XM*RTEMP2) 

MAIARY(IDXZY,J,I) = (YM*RTEMP2) 

MAIARY ( IDXMSK, J , I ) = 1. 

450 WRITE(18, 20520) MAIARY (IDXZ , J, I) , MAIARY (IDXZX ,J,I) 

1 , MAIARY (IDXZY ,J, I) ,MAIARY( IDXMSK , J, I) 

20800 CONTINUE 
20900 CONTINUE 
END 


DRS04610 

DRS04620 

=DRS04630 

DRS04640 

DRS04650 

DRS04660 

DRS04670 

DRS04680 

DRS04690 

DRS04700 

DRS04710 

DRS04720 

DRS04730 

DRS04740 
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PROGRAM DUALREF DUAOOOIO 

********* *********************************************************** DUA00020 

* AUTHOR : R. ACOSTA AND A. LAGIN VM VERSION DUA00030 

* DATE : 7/15/91 DUA00040 

* PURPOSE : TO COMPUTE FAR FIELD CO-POL AND CROSS POL OF A DUA00050 

* GENERALIZED DUAL REFLECTOR SYSTEM DUA00060 

* GENERALIZED DUAL REFLECTOR ANALYSIS PROGRAM DUA00070 


****************************************************************** ******Q(J^QQ080 

** COMPILE CONSTANTS DUA00090 


RFAI f 

PARAMETER (C 
REAL PI 

PARAMETER (PI 
REAL ETA 

PARAMETER (ETA 
REAL R2DEG 

PARAMETER (R2DEG 
REAL D2RAD 

PARAMETER (D2RAD 
INTEGER IDXZ 

INTEGER IDXZX 

INTEGER IDXZY 

INTEGER IDXMSK 

INTEGER IDXVCX 

INTEGER IDXVCY 

INTEGER IDXVCZ 

INTEGER IDXRMJ 

INTEGER IDXIMJ 

INTEGER IDXUNM 

INTEGER IDXJVX 

INTEGER IDXNRM 

INTEGER IDXDTX 

INTEGER IDXMNM 

INTEGER IDXMSI 

INTEGER IDXAOT 

INTEGER IDXPWR 

INTEGER IDXADB 

INTEGER IDXRDB 

INTEGER MXMX , MXMY , TMXMX , TMXMY 

INTEGER MXSX,MXSY,TMXSX,TMXSY 

INTEGER MXTHE.MXPHI 

INTEGER MXYZ 

INTEGER DXYZ 

INTEGER LNRY 

INTEGER VDTX 

INTEGER NDTX 

INTEGER MDTX 

INTEGER KDTX 

PARAMETER (IDXZ = 1 
PARAMETER (IDXZX = 2 
PARAMETER (IDXZY = 3 
PARAMETER (IDXMSK = 4 
PARAMETER (IDXVCX = 1 
PARAMETER (IDXVCY = 2 
PARAMETER (IDXVCZ = 3 
PARAMETER (IDXRMJ = 1 
PARAMETER (IDXIMJ = 2 
PARAMETER (IDXUNM = 3 
PARAMETER (IDXJVX = 1 


DUAOOllO 
DUA00120 
DUA00130 
DUA00140 
DUA00150 
DUA00160 
DUA00170 
DUA00180 
DUA00190 
DUA00200 
DUA00210 
DUA00220 
DUA00230 
DUA00240 
DUA00250 
DUA00260 
DUA00270 
DUA00280 
DUA00290 
DUA00300 
DUA00310 
DUA00320 
DUA00330 
DUA00340 
DUA00350 
DUA00360 
DUA00370 
DUA00380 
DUA00390 
DUA00400 
DUA00410 
DUA00420 
DUA00430 
DUA00440 
DUA00450 
DUA00460 
DUA00470 
DUA00480 
DUA00490 

) DUA00500 
) DUA00510 
) DUA00520 
) DUA00530 
) DUA00540 
) DUA00550 
) DUA00560 
) DUA00570 
) DUA00580 
) DUA00590 
) DUA00600 


= 2.99792458E+08 ) 
= 3.141592653589793238) 
= 4.0*PI*C*(1 .E-7) ) 
= 180. /PI ) 

= PI/180. ) 
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PARAMETER (IDXNRM = 2 
PARAMETER (IDXDTX = 3 
PARAMETER (IDXMNM = 1 
PARAMETER (IDXMSI = 2 
PARAMETER (IDXAOT = 3 
PARAMETER (IDXPWR = 1 
PARAMETER (IDXADB = 2 

PARAMETER (IDXRDB = 3 ) 

********** ************************************* ******************** 

***** jNPUT : NUMBER OF POINTS IN X AND Y MAIN REF. GRID.*********** 
******************************************************************* 


PARAMETER (MXMX = 101 ) 

PARAMETER (MXMY = 101 ) 

********************************************************************** 

****** INPUT . NUMBER OF POINTS IN X AND Y SUB REF. GRID*************** 
********************************************************************** 

PARAMETER (MXSX * 61 ) 

PARAMETER (MXSY = 61 ) 

********************************************************************** 


) DUA00610 

) DUA00620 

) DUA00630 

) DUA00640 

DUA00650 
DUA00660 
DUA00670 
DUA00680 
DUA00690 
DUA00700 
DUA00710 
DUA00720 
DUA00730 
DUA00740 
DUA00750 
DUA00760 
DUA00770 
DUA00780 
DUA00790 

) DUA00800 

) DUA00810 

) DUA00820 

) DUA00830 

DUA00840 
DUA00850 
DUA00860 
DUA00870 
DUA00880 
DUA00890 

) DUA00900 

) DUA00910 

) DUA00920 

) DUA00930 

) DUA00940 

) DUA00950 

) DUA00960 

************************************************************************qu/\qo970 
** INTRINSIC FUNCTIONS DUA00980 

********************************** **************************************QJJ^QQQQQ 

INTRINSIC SQRT DUA01000 

INTRINSIC SIN DUA01010 

INTRINSIC COS DUA01020 

INTRINSIC ACOS DUA01030 

INTRINSIC NINT DUA01040 

************************************************************************0U/\giQ50 
** EXTERNAL FUNCTIONS DUA01060 


PARAMETER (TMXMX 
PARAMETER (TMXMY 
PARAMETER (TMXSX 
PARAMETER (TMXSY 


= MXMX + 1 
= MXMY + 1 
= MXSX + 1 
= MXSY + 1 


********************************************************************** 

******j^pyj . number OF FAR-FIELD grid points************************* 
********************************************************************** 

PARAMETER (MXTHE = 50 ) 

PARAMETER (MXPHI = 360 ) 

********************************************************************** 


PARAMETER (MXYZ 

= 4 

PARAMETER (DXYZ 

= 3 

PARAMETER (LNRY 

= 1 

PARAMETER (VDTX 

= 4 

PARAMETER (NDTX 

= 3 

PARAMETER (MDTX 

= 10 

PARAMETER (KDTX 

= 18 


******************************************** ****************************0(J^Q J Q7Q 


REAL DOT DUA01080 

EXTERNAL DOT DUA01090 

REAL FDPTRN DUA01100 

EXTERNAL FDPTRN DUA01110 

************************************************************************00^01 120 

** EXTERNAL SUBROUTINES DUA01130 

******************************* ************************************** ***ou/\o 1140 
EXTERNAL CROSS DUA01150 

EXTERNAL SCALER DUA01160 

EXTERNAL VECADD DUA01170 

EXTERNAL VECSUB DUA01180 

************************************************************************00^0 ^ y 90 
** RUN TIME CONSTANTS DUA01200 
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************************************************************************^^^210 
REAL DTAARY (MDTX) DUA01220 


REAL 

REAL 

REAL 

REAL 

REAL 

INTEGER 

INTEGER 

REAL 

REAL 


DTAARY (MDTX) 
DMX.DMY 
INCMX, INCMY 
DSX.DSY 
INCSX, INCSY 
MAXMX, MAXMY 
MAXSX, MAXSY 
FREQ 
Q 


DUA01230 

DUA01240 

DUA01250 

DUA01260 

DUA01270 

DUA01280 

DUA01290 

DUA01300 


************* ****** ****** ***********************************************QU^Q 13 10 

** EQUIVALENCE DUA01320 

************************************************************************ 0 U ^ O 1230 

EQUIVALENCE (DTAARY ( 1),DMX ) DUA01340 

EQUIVALENCE (DTAARY ( 2) t DMY ) DUA01350 

EQUIVALENCE (DTAARY ( 3), INCMX) DUA01360 

EQUIVALENCE (DTAARY ( 4), INCMY) DUA01370 

EQUIVALENCE (DTAARY ( 5),DSX ) DUA01380 

EQUIVALENCE (DTAARY ( 6),DSY ) DUA01390 

EQUIVALENCE (DTAARY ( 7), INCSX) DUA01400 

EQUIVALENCE (DTAARY ( 8) , INCSY) DUA01410 

EQUIVALENCE (DTAARY ( 9), FREQ ) DUA01420 

EQUIVALENCE (DTAARY(IO) ,Q ) DUA01430 

************************************************************************00^01440 

** RUN TIME VARIABLES DUA01450 

********************************************************************** **qUAO 1460 

REAL LAMBDA DUA01470 

REAL K DUA01480 

*********************************************************************** *pu/\oi490 
** EQUIVALENCE DUA01500 

********************** **************************************************Q\jf{Qi§iQ 

** Input Arrays DUA01520 

** MAIFIL ==> MAIARY() DUA01530 

** SUBFIL ==> SUBARY () DUA01540 

** XYZFIL ==> XYZARY() ==> MRXYZO() ,SRXYZO() f FEDXYZ() ,REFXYZ() DUA01550 

** DTAFIL ==> DTAARYQ ==> DMX.DMY, INCMX, INCMY.DSX.DSY, INCSX, INCSY, FREQDUA01560 

********************** **************************************************Q\j^Q^ 570 

******READ IN THE XYZARY AND DTAARY FROM FILE GENERATOR*****************DUA01580 
******************************* *** *** ********************************** *QUA0 1 590 


DO 10 1=1,10 
READ(15,805) DTAARY(I) 

805 FORMAT (5X, E15.8) 

10 CONTINUE 

***************************************************************** 
******************************************************************* 
** Initial Calculations 

****************************************************************** 
LAMBDA = C/FREQ 
K = 2*PI/LAMBDA 
MAXMX = (NINT (DMX/INCMX) + 1) 

MAXMY = (NINT (DMY/INCMY) + 1) 

MAXSX = (NINT (DSX/INCSX) + 1) 

MAXSY = (NINT (DSY/INCSY) + 1) 

CALL SUBMAI (MAXMX, MAXMY, MAXSX, MAXSY, MDTX, DTAARY) 

END 

******************************************************************* 

************^ PROGRAM******************************************* 

******************************************************************* 

SUBROUTINE SUBMAI ( PMXMX , PMXMY , PMXSX , PMXSY, PDTX , DTXARX) 

INTEGER PMXMX 


DUA01600 

DUA01610 

DUA01620 

DUA01630 

DUA01640 

DUA01650 

DUA01660 

DUA01670 
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INTEGER PMXMY 

INTEGER PMXSX 

INTEGER PMXSY 

INTEGER PDTX 

REAL DTXARX(PDTX) 

**************************************************** 


** COMPILE CONSTANTS 

**************************************************** 


REAL 

PARAMETER 

REAL 

PARAMETER 

REAL 

PARAMETER 

REAL 

PARAMETER 

REAL 

PARAMETER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

INTEGER 


C 

(C 

PI 

(PI 

ETA 

(ETA 

R2DEG 

(R2DEG 

D2RAD 

(D2RAD 

IDXZ 

IDXZX 

IDXZY 

IDXMSK 

IDXVCX 

IDXVCY 

IDXVCZ 

IDXRMJ 

IDXIMJ 

IDXUNM 

IDXJVX 

IDXNRM 

IDXDTX 

IDXMNM 

IDXMSI 

IDXAOT 

IDXPWR 


= 2.99792458E+08 


= 3.141592653589793238) 
= 4.0*PI*C*(l.E-7) ) 

= 180. /PI ) 

= PI/180. ) 


DUA01810 

DUA01820 

DUA01830 

DUA01840 

DUA01850 

********************Qy^Q IggQ 

DUA01870 

***************** ***[)u/\q j 880 

DUA01890 
DUA01900 
DUA01910 
DUA01920 
DUA01930 
DUA01940 
DUA01950 
DUA01960 
DUA01970 
DUA01980 
DUA01990 
DUA02000 
DUA02010 
DUA02020 
DUA02030 
DUA02040 
DUA02050 
DUA02060 
DUA02070 
DUA02080 
DUA02090 
DUA02100 
DUA021 10 
DUA02120 
DUA02130 
DUA02140 
DUA02150 


INTEGER IDXADB 

INTEGER IDXRDB 

INTEGER MXMX , MXMY , TMXMX , TMXMY 

INTEGER MXSX,MXSY f TMXSX,TMXSY 

INTEGER MXTHE.MXPHI 

INTEGER MXYZ 

INTEGER DXYZ 

INTEGER LNRY 

INTEGER VDTX 

INTEGER NDTX 

INTEGER MDTX 

INTEGER KDTX 


PARAMETER (IDXZ = 1 
PARAMETER (IDXZX = 2 
PARAMETER (IDXZY = 3 
PARAMETER (IDXMSK = 4 
PARAMETER (IDXVCX = 1 
PARAMETER (IDXVCY = 2 
PARAMETER (IDXVCZ = 3 
PARAMETER (IDXRMJ = 1 
PARAMETER (IDXIMJ = 2 
PARAMETER (IDXUNM = 3 
PARAMETER (IDXJVX = 1 
PARAMETER (IDXNRM = 2 
PARAMETER (IDXDTX = 3 


) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 


DUA02160 

DUA02170 

DUA02180 

DUA02190 

DUA02200 

DUA02210 

DUA02220 

DUA02230 

DUA02240 

DUA02250 

DUA02260 

DUA02270 

DUA02280 

DUA02290 

DUA02300 

DUA02310 

DUA02320 

DUA02330 

DUA02340 

DUA02350 

DUA02360 

DUA02370 

DUA02380 

DUA02390 

DUA02400 
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PARAMETER (IDXMNM = 1 ) 

PARAMETER (IDXMSI = 2 ) 

PARAMETER (IDXAOT = 3 ) 

PARAMETER (IDXPWR = 1 ) 

PARAMETER (IDXADB = 2 ) 

PARAMETER (IDXRDB = 3 ) 

********************************************************************* 

****** INPUT . NUMBER OF POINTS IN X AND Y MAIN REF. GRID************* 
********************************************************************* 

PARAMETER (MXMX = 101 ) 

PARAMETER (MXMY = 101 ) 

********************************************************************* 

****** INPUT . NUMBER OF POINTS IN X AND Y SUB REF. GRID************** 
********************************************************************* 

PARAMETER (MXSX = 61 ) 

PARAMETER (MXSY = 61 ) 

********************************************************************* 

PARAMETER (TMXMX = MXMX + 1 ) 

= MXMY + 1 ) 

= MXSX + 1 ) 

= MXSY +1 ) 

******************************************************************** 

******* INPUT : FAR-FIELD GRID POINTS THETA AND pm***************** 
■***■*■******** *********** ********** *********★*******'******************** 

PARAMETER (MXTHE = 50 ) 

PARAMETER (MXPHI = 360 ) 

******** ****************** ***************************************** 


PARAMETER (TMXMY 
PARAMETER (TMXSX 
PARAMETER (TMXSY 


DUA02410 

DUA02420 

DUA02430 

DUA02440 

DUA02450 

DUA02460 

DUA02470 

DUA02480 

DUA02490 

DUA02500 

DUA02510 

DUA02520 

DUA02530 

DUA02540 

DUA02550 

DUA02560 

DUA02570 

DUA02580 

DUA02590 

DUA02600 

DUA02610 

DUA02620 

DUA02630 

DUA02640 

DUA02650 

DUA02660 

DUA02670 


PARAMETER (MXYZ 

= 4 

) 

DUA02680 

PARAMETER (DXYZ 

= 3 

) 

DUA02690 

PARAMETER (LNRY 

= 1 

) 

DUA02700 

PARAMETER (VDTX 

= 4 

) 

DUA02710 

PARAMETER (NDTX 

= 3 

) 

DUA02720 

PARAMETER (MDTX 

= 10 

) 

DUA02730 

PARAMETER (KDTX 

= 18 

) 

DUA02740 


************************************************************************QUA02750 
** INTRINSIC FUNCTIONS DUA02760 

************************************************************************0UAO277O 
INTRINSIC SQRT DUA02780 

INTRINSIC SIN DUA02790 

INTRINSIC COS DUA02800 

INTRINSIC ACOS DUA02810 

INTRINSIC NINT DUA02820 

************************************************************************0UAO283O 
** EXTERNAL FUNCTIONS DUA02840 

******************************************************************** ****0UAO285O 
REAL DOT DUA02860 

EXTERNAL DOT DUA02870 

REAL FDPTRN DUA02880 

EXTERNAL FDPTRN DUA02890 

*********************************************************** *************0UAO29OO 
** EXTERNAL SUBROUTINES DUA02910 

***************★★★*★**********************★**★***★*★*★*★***★**★*********[) U AO 2920 
EXTERNAL CROSS DUA02930 

EXTERNAL SCALER DUA02940 

EXTERNAL VECADD DUA02950 

EXTERNAL VECSUB DUA02960 

********************************************* ***************************0U^O29 70 

** RUN TIME CONSTANTS DUA02980 

****************************************************************** ******0UAO299O 
REAL MAIARY (VDTX, TMXMY, TMXMX) DUA03000 
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REAL SUBARY (VDTX, TMXSY, TMXSX) 

REAL XYZARY (DXYZ.MXYZ) 

REAL SRXYZO(DXYZ) ,XSO,YSO,ZSO 

REAL MRXYZO(DXYZ) ,XMO,YMO r ZMO 

REAL FEDXYZ(DXYZ) ,XF , YF ,ZF 

REAL REFXYZ (DXYZ) , XR , YR ,ZR 

REAL SR(DXYZ) , SRX ,SRY ,SRZ 

REAL USR(DXYZ) , USRX, USRY, USRZ 

REAL MAGSR 

REAL DTAARY (MDTX) 

REAL DMX.DMY 

REAL INCMX, INCMY 

REAL DSX.DSY 

REAL INCSX, INCSY 

REAL FREQ 

REAL Q 

REAL DTXARY(KDTX) 

REAL BEGPHX 

REAL ENDPHX 

REAL IDXPHX 

REAL STPPHX 

REAL INCPHX 

REAL BEGTHX 

REAL ENDTHX 

REAL IDXTHX 

REAL STPTHX 

REAL INCTHX 

REAL PTTMNX 

REAL PTTMXX 

REAL ADBMNX 

REAL ADBMXX 

REAL RDBMNX 

REAL PRAD 

REAL RINTNS 

REAL RIFCTR 

REAL DIRCTV 

********************************************************************* 

** EQUIVALENCE 

******************************************************************** 
EQUIVALENCE (XYZARY(1 , 1) ,MRXYZO(l) ) 

EQUIVALENCE (XYZARY(1 , 2) , SRXYZO(l) ) 

EQUIVALENCE (XYZARY(1 , 3) , FEDXYZ (1) ) 

EQUIVALENCE (XYZARY (1 , 4) , REFXYZ (1) ) 

EQUIVALENCE (MRXYZO( l),XMO ) 

EQUIVALENCE (MRXYZO ( 2),YM0 ) 

EQUIVALENCE (MRXYZO( 3),ZMO ) 

EQUIVALENCE (SRXYZO( l),XSO ) 

EQUIVALENCE (SRXYZO( 2),YS0 ) 

EQUIVALENCE (SRXYZO( 3),ZSO ) 

EQUIVALENCE (FEDXYZ ( 1),XF ) 

EQUIVALENCE ( FEDXYZ ( 2),YF ) 

EQUIVALENCE (FEDXYZ ( 3),ZF ) 

EQUIVALENCE (REFXYZ ( 1) , XR ) 

EQUIVALENCE (REFXYZ ( 2),YR ) 

EQUIVALENCE (REFXYZ ( 3),ZR ) 

EQUIVALENCE (SR ( 1) , SRX ) 

EQUIVALENCE (SR ( 2),SRY ) 

EQUIVALENCE (SR ( 3),SRZ ) 

EQUIVALENCE (USR ( 1),USRX ) 

EQUIVALENCE (USR ( 2), USRY ) 


DUA03010 

DUA03020 

DUA03030 

DUA03040 

DUA03050 

DUA03060 

DUA03070 

DUA03080 

DUA03090 

DUA03100 

DUA031 10 

DUA03120 

DUA03130 

DUA03140 

DUA03150 

DUA03160 

DUA03170 

DUA03180 

DUA03190 

DUA03200 

DUA03210 

DUA03220 

DUA03230 

DUA03240 

DUA03250 

DUA03260 

DUA03270 

DUA03280 

DUA03290 

DUA03300 

DUA03310 

DUA03320 

DUA03330 

DUA03340 

DUA03350 

DUA03360 

DUA03370 

DUA03380 

DUA03390 

DUA03400 

DUA03410 

DUA03420 

DUA03430 

DUA03440 

DUA03450 

DUA03460 

DUA03470 

DUA03480 

DUA03490 

DUA03500 

DUA03510 

DUA03520 

DUA03530 

DUA03540 

DUA03550 

DUA03560 

DUA03570 

DUA03580 

DUA03590 

DUA03600 
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EQUIVALENCE 

(USR ( 3) , USRZ ) 

DUA03610 

EQUIVALENCE 

(DTAARY ( 1 ) , DMX ) 

DUA03620 

EQUIVALENCE 

(DTAARY ( 2) , DMY ) 

DUA03630 

EQUIVALENCE 

(DTAARY ( 3) , INCMX) 

DUA03640 

EQUIVALENCE 

(DTAARY ( 4) , INCMY) 

DUA03650 

EQUIVALENCE 

(DTAARY ( 5) , DSX ) 

DUA03660 

EQUIVALENCE 

(DTAARY ( 6) , DSY ) 

DUA03670 

EQUIVALENCE 

(DTAARY ( 7) , INCSX) 

DUA03680 

EQUIVALENCE 

(DTAARY ( 8) , INCSY) 

DUA03690 

EQUIVALENCE 

(DTAARY ( 9), FREQ ) 

DUA03700 

EQUIVALENCE 

(DTAARY(IO) ,Q ) 

DUA03710 

EQUIVALENCE 

(DTXARY ( 1) , BEGPHX) 

DUA03720 

EQUIVALENCE 

(DTXARY ( 2) , ENDPHX) 

DUA03730 

EQUIVALENCE 

(DTXARY ( 3) , IDXPHX) 

DUA03740 

EQUIVALENCE 

(DTXARY ( 4) , STPPHX) 

DUA03750 

EQUIVALENCE 

(DTXARY ( 5) , INCPHX) 

DUA03760 

EQUIVALENCE 

(DTXARY ( 6) , BEGTHX) 

DUA03770 

EQUIVALENCE 

(DTXARY ( 7) , ENDTHX) 

DUA03780 

EQUIVALENCE 

(DTXARY ( 8) , IDXTHX) 

DUA03790 

EQUIVALENCE 

(DTXARY ( 9) , STPTHX) 

DUA03800 

EQUIVALENCE 

(DTXARY (10) , INCTHX) 

DUA03810 

EQUIVALENCE 

(DTXARY(ll) , PTTMNX) 

DUA03820 

EQUIVALENCE 

(DTXARY (12) ,PTTMXX) 

DUA03830 

EQUIVALENCE 

(DTXARY(13) , ADBMNX) 

DUA03840 

EQUIVALENCE 

(DTXARY ( 14) r ADBMXX) 

DUA03850 

EQUIVALENCE 

(DTXARY(15) f RDBMNX) 

DUA03860 

EQUIVALENCE 

(DTXARY(16) , PRAD ) 

DUA03870 

EQUIVALENCE 

(DTXARY(17) ,RINTNS) 

DUA03880 

EQUIVALENCE 

(DTXARY(18) ,DIRCTV) 

DUA03890 


******************************** ****************************************qu^q 3 goo 
** RUN TIME VARIABLES DUA03910 
************************ *********** *************************************qu ao 3 920 


REAL 

MCDARY (DXYZ , DXYZ , TMXMY , TMXMX) 

DUA03930 

REAL 

SCDARY (DXYZ , DXYZ, TMXSY , TMXSX) 

DUA03940 

REAL 

MAIXYZ (DXYZ) , XM, YM, ZM 

DUA03950 

REAL 

SUBXYZ(DXYZ) ,XS, YS, ZS 

DUA03960 

REAL 

GENXYZ(DXYZ) , X I , YJ, ZIJ 

DUA03970 

REAL 

TMPXYZ(DXYZ,2) 

DUA03980 

REAL 

TMRXYZ(DXYZ) , TMRX, TMRY, TMRZ 

DUA03990 

REAL 

TMIXYZ (DXYZ) , TMIX, TMIY, TMIZ 

DUA04000 

REAL 

HFLD(DXYZ,2) 

DUA04010 

REAL 

HVR (DXYZ) ,HVRX,HVRY,HVRZ 

DUA04020 

REAL 

HVI (DXYZ) ,HVIX,HVIY,HVIZ 

DUA04030 

REAL 

JFLD(DXYZ,2) 

DUA04040 

REAL 

JVR (DXYZ) , JVRX, JVRY, JVRZ 

DUA04050 

REAL 

JVI (DXYZ),JVIX,JVIY,JVIZ 

DUA04060 

REAL 

SUM (DXYZ, 2) 

DUA04070 

REAL 

SUMR(DXYZ) ,SUMRX,SUMRY,SUMRZ 

DUA04080 

REAL 

SUMI (DXYZ) ,SUMIX,SUMIY,SUMIZ 

DUA04090 

REAL 

NORM(DXYZ) , NX , NY ,NZ 

DUA04100 

REAL 

MAGNRM 

DUA04110 

REAL 

SI (DXYZ), SIX , SIY , SIZ 

DUA04120 

REAL 

USI (DXYZ) , US IX, US I Y, US I Z 

DUA04130 

REAL 

MAGS I 

DUA04140 

REAL 

PV (DXYZ) , PVX , PVY ,PVZ 
UPV (DXYZ) ,UPVX,UPVY,UPVZ 

DUA04150 

REAL 

DUA04160 

REAL 

MAGPV 

DUA04170 

REAL 

HV (DXYZ) , HVX ,HVY ,HVZ 

DUA04180 

REAL 

UHV (DXYZ) ,UHVX,UHVY,UHVZ 

DUA04190 

REAL 

MAGHV 

DUA04200 
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REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

INTEGER 

REAL 

REAL 

REAL 

INTEGER 

INTEGER 

REAL 

REAL 

REAL 

INTEGER 

REAL 

REAL 

REAL 

INTEGER 

INTEGER 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

REAL 

INTEGER 

INTEGER 

REAL 

REAL 

REAL 

INTEGER 

INTEGER 

INTEGER 

INTEGER 

CHARACTER 

REAL AAAKR 


JV (DXYZ) , JVX , JVY , JVZ 
UJV(DXYZ) ,UJVX,l)JVY,UJVZ 
MAGJV 

R1 (DXYZ).RIX ,R1Y ,R1Z 
UR1 (DXYZ) ,URlX t URlY,URlZ 
MAGR1 

INTG(DXYZ) , INTX, INTY, INTZ 
MAGINT 

RFF (DXYZ) , RFFX, RFFY, RFFZ 

PTTRN (NDTX , 0 : MXTHE, 0 :MXPHI) 

PTTMIN 

PTTMAX 

LAMBDA 

K 

KR 

RR 

PSI (2) , COSKR, SINKR 

MIN 

MAX 

R1TMP 

R2TMP 

CMPTMP(2) , CMPTMR.CMPTMI 

SCALE 

ANGLE 

ANGPHX 

ANGPHI 

BEGPHI 

ENDPHI 

IDXPHI 

STPPHI 

INCPHI 

SINPHI 

COSPHI 

ANGTHX 

ANGTHE 

BEGTHE 

ENDTHE 

IDXTHE 

STPTHE 

INCTHE 

SINTHE 

COSTHE 

EPHRE 

EPHIM 

ETHRE 

ETHIM 

COZ 

MSKFAC 
MAXMX , MAXMY 
MAXSX.MAXSY 
DS 

INCMXY 

INCSXY 

I,J,IP,JP,V f W 

IOS 

FLG 

ITMP 

TIME*8 


INTEGER III KR 


DUA04210 

DUA04220 

DUA04230 

DUA04240 

DUA04250 

DUA04260 

DUA04270 

DUA04280 

DUA04290 

DUA04300 

DUA04310 

DUA04320 

DUA04330 

DUA04340 

DUA04350 

DUA04360 

DUA04370 

DUA04380 

DUA04390 

DUA04400 

DUA04410 

DUA04420 

DUA04430 

DUA04440 

DUA04450 

DUA04460 

DUA04470 

DUA04480 

DUA04490 

DUA04500 

DUA04510 

DUA04520 

DUA04530 

DUA04540 

DUA04550 

DUA04560 

DUA04570 

DUA04580 

DUA04590 

DUA04600 

DUA04610 

DUA04620 

DUA04630 

DUA04640 

DUA04650 

DUA04660 

DUA04670 

DUA04680 

DUA04690 

DUA04700 

DUA04710 

DUA04720 

DUA04730 

DUA04740 

DUA04750 

DUA04760 

DUA04770 

DUA04780 

DUA04790 

DUA04800 
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t***********************************************************************^^^^ 

** EQUIVALENCE DUA04820 

***************** t******************************************************^^^#^ 


EQUIVALENCE 

(SUBXYZ(l) ,XS ) 

DUA04840 

EQUIVALENCE 

(SUBXYZ(2) , YS ) 

DUA04850 

EQUIVALENCE 

(SUBXYZ (3) , ZS ) 

DUA04860 

EQUIVALENCE 

(MAIXYZ(l) , XM ) 

DUA04870 

EQUIVALENCE 

(MAIXYZ(2) , YM ) 

DUA04880 

EQUIVALENCE 

(MAIXYZ(3) ,ZM ) 

DUA04890 

EQUIVALENCE 

(GENXYZ(l) ,XI ) 

DUA04900 

EQUIVALENCE 

(GENXYZ(2) , YJ ) 

DUA04910 

EQUIVALENCE 

(GENXYZ (3) , ZIJ ) 

DUA04920 

EQUIVALENCE 

(TMPXYZ(l , 1) ,TMRXYZ(l) ) 

DUA04930 

EQUIVALENCE 

(TMPXYZ(1,2) , TMIXYZ ( 1 ) ) 

DUA04940 

EQUIVALENCE 

(TMRXYZ(l) , TMRX ) 

DUA04950 

EQUIVALENCE 

(TMRXYZ(2) , TMRY ) 

DUA04960 

EQUIVALENCE 

(TMRXYZ(3) ,TMRZ ) 

DUA04970 

EQUIVALENCE 

(TMIXYZ (1) , TMIX ) 

DUA04980 

EQUIVALENCE 

(TMIXYZ (2) , TMIY ) 

DUA04990 

EQUIVALENCE 

(TMIXYZ (3), TMIZ ) 

DUA05000 

EQUIVALENCE 

(SUM (1,1) , SUMR ( 1 ) ) 

DUA05010 

EQUIVALENCE 

(SUM(1 , 2) , SUMI (1) ) 

DUA05020 

EQUIVALENCE 

(SUMR 

(1) ,SUMRX ) 

DUA05030 

EQUIVALENCE 

(SUMR 

(2) , SUMRY ) 

DUA05040 

EQUIVALENCE 

(SUMR 

(3) , SUMRZ ) 

DUA05050 

EQUIVALENCE 

(SUMI 

(l).SUMIX ) 

DUA05060 

EQUIVALENCE 

(SUMI 

(2) , SUMIY ) 

DUA05070 

EQUIVALENCE 

(SUMI 

(3) ,SUMIZ ) 

DUA05080 

EQUIVALENCE 

(NORM 

(1) » NX ) 

DUA05090 

EQUIVALENCE 

(NORM 

(2), NY ) 

DUA05100 

EQUIVALENCE 

(NORM 

(3) ,NZ ) 

DUA05110 

EQUIVALENCE 

(SI 

(1) t SIX ) 

DUA05120 

EQUIVALENCE 

(SI 

(2) , SIY ) 

DUA05130 

EQUIVALENCE 

(SI 

(3) ,SIZ ) 

DUA05140 

EQUIVALENCE 

(USI 

(1) rUSIX ) 

DUA05150 

EQUIVALENCE 

(USI 

(2) ,USIY ) 

DUA05160 

EQUIVALENCE 

(USI 

(3) , USIZ ) 

DUA05170 

EQUIVALENCE 

(PV 

(D,PVX ) 

DUA05180 

EQUIVALENCE 

(PV 

(2) , PVY ) 

DUA05190 

EQUIVALENCE 

(PV 

(3) , PVZ ) 

DUA05200 

EQUIVALENCE 

(UPV 

(1),UPVX ) 

DUA05210 

EQUIVALENCE 

(UPV 

(2) ,UPVY ) 

DUA05220 

EQUIVALENCE 

UPV 

(3) . UPVZ ) 

DUA05230 

EQUIVALENCE 

(HFLD(1,1),HVR(1)) 

DUA05240 

EQUIVALENCE 

(HFLD(1 , 2) ,HVI (1) ) 

DUA05250 

EQUIVALENCE 

(HVR 

(1),HVRX ) 

DUA05260 

EQUIVALENCE 

(HVR 

(2) , HVRY ) 

DUA05270 

EQUIVALENCE 

(HVR 

(3) f HVRZ ) 

DUA05280 

EQUIVALENCE 

(HVI 

(1) ,HVIX ) 

DUA05290 

EQUIVALENCE 

(HVI 

(2) , HVIY ) 

DUA05300 

EQUIVALENCE 

(HVI 

(3) ,HVIZ ) 

DUA05310 

EQUIVALENCE 

(HV 

(1) , HVX ) 

DUA05320 

EQUIVALENCE 

(HV 

(2) , HVY ) 

DUA05330 

EQUIVALENCE 

(HV 

(3) , HVZ ) 

DUA05340 

EQUIVALENCE 

(UHV 

(l).UHVX ) 

DUA05350 

EQUIVALENCE 

(UHV 

(2) ,UHVY ) 

DUA05360 

EQUIVALENCE 

(UHV 

(3) ,UHVZ ) 

DUA05370 

EQUIVALENCE 

( J FLD ( 1 , 1) , JVR(l) ) 

DUA05380 

EQUIVALENCE 

( J FLD (1,2) , JVI(l)) 

DUA05390 

EQUIVALENCE 

(JVR 

(1) , JVRX ) 

DUA05400 
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EQUIVALENCE 

(JVR 

(2) , JVRY ) 

DUA05410 

EQUIVALENCE (JVR 

(3) , JVRZ ) 

DUA05420 

EQUIVALENCE 

(JVI 

(i),ovix ) 

DUA05430 

EQUIVALENCE 

(JVI 

(2) , JVI Y ) 

DUA05440 

EQUIVALENCE 

(JVI 

(3) , JVIZ ) 

DUA05450 

EQUIVALENCE 

(J V 

(1) , JVX ) 

DUA05460 

EQUIVALENCE 

(JV 

(2) , JVY ) 

DUA05470 

EQUIVALENCE 

(JV 

(3) , JVZ ) 

DUA05480 

EQUIVALENCE 

(UJV 

(1),UJVX ) 

DUA05490 

EQUIVALENCE 

(UJV 

(2) , UJVY ) 

DUA05500 

EQUIVALENCE 

(UJV 

(3) , UJVZ ) 

DUA05510 

EQUIVALENCE 

(R1 

(1),R1X ) 

DUA05520 

EQUIVALENCE 

(R1 

(2) , R1Y ) 

DUA05530 

EQUIVALENCE 

(R1 

(3) , R1Z ) 

DUA05540 

EQUIVALENCE 

( UR 1 

(1),UR1X ) 

DUA05550 

EQUIVALENCE 

(UR1 

(2) , UR1 Y ) 

DUA05560 

EQUIVALENCE 

(UR1 

(3) , UR1Z ) 

DUA05570 

EQUIVALENCE 

(INTG 

U),INTX ) 

DUA05580 

EQUIVALENCE 

(INTG 

(2) , INTY ) 

DUA05590 

EQUIVALENCE 

(INTG 

(3) , INTZ ) 

DUA05600 

EQUIVALENCE 

(RFF 

(1) , RFFX ) 

DUA05610 

EQUIVALENCE 

(RFF 

(2) , RFFY ) 

DUA05620 

EQUIVALENCE 

(RFF 

(3) , RFFZ ) 

DUA05630 

EQUIVALENCE 

(PSI 

(1) ,COSKR) 

DUA05640 

EQUIVALENCE 

(PSI 

(2) , SINKR) 

DUA05650 

EQUIVALENCE 

(CMPTMP(l) ,CMPTMR) 

DUA05660 

EQUIVALENCE 

(CMPTMP(2) .CMPTMI) 

DUA05670 


* ********************************************************************** *DUAO 5680 
** INITIALIZE INPUT ARRAYS DUA05690 

** MAIFIL ==> MAIARY () DUA05700 

** SUBFIL ==> SUBARY () DUA05710 

** XYZFIL ==> XYZARY () ==> MRXYZOQ ,SRXYZO() ,FEDXYZ() ,REFXYZ() DUA05720 

** DTAFIL ==> DTAARY() ==> DMX.DMY, INCMX, INCMY.DSX.DSY, INCSX, INCSY, FREQ,DUA05730 

************************************************************************pjj^Q^-^Q 

** GENERATE OUTPUT ARRRAYS DUA05750 

** MCDFIL <== MCDARY () DUA05760 

** SCDFIL <= = SCDARY () DUA05770 

** RFXFIL <== PTTRN () DUA05780 

** OTXFIL <== DTXARY () <== BEG, END, IDX.STP, INC * THE, PHI , PTT * ADR * DUA05790 

********************************************** **************************py^Qj-gQQ 

******rEa D in the XYZARY, SUBARY, MAIARY ARRAYS FROM FILE GEN.************DUA05810 
************************************************************************^ UA05820 


DO 11 11=1,4 DUA05830 

DO 12 JJ=1 , 3 DUA05840 

READ (16, 806) XYZARY ( J J ,11) DUA05850 

806 FORMAT (5X, E15.8) DUA05860 

12 CONTINUE DUA05870 

11 CONTINUE DUA05880 


***************************** *************************************** ****Q[J^Q£jggQ 


***************************************************************** 
DO 13 I I I = 1 , TMXSX, 1 
DO 14 J J J = 1 , TMXSY, 1 
READ (17, 807 ) SUBARY ( IDXZ , J J J ,111) , 

1 SUBARY (IDXZX, JJJ, III) , 

2 SUBARY ( IDXZ Y , JJJ , I I I ) , 

3 SUBARY (IDXMSK, JJJ, III) 

807 FORMAT (5X,4(E15.8,2X)) 

14 CONTINUE 

13 CONTINUE 

DO 15 1 1 1 1 = 1 , TMXMX, 1 


DUA05900 

DUA05910 

DUA05920 

DUA05930 

DUA05940 

DUA05950 

DUA05960 

DUA05970 

DUA05980 

DUA05990 

DUA06000 
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DO 16 JJJJ=1,TMXMY,1 DUA06010 

READ ( 18 1 807 ) MAIARY ( IDXZ , J J J J ,1111), DUA06020 

1 MAIARY(IDXZX, JJJJ, IIII) , DUA06030 

2 MAIARY(IDXZY, JJJJ, IIII), DUA06040 

3 MAIARY(IDXMSK, JJJJ, IIII) DUA06050 

16 CONTINUE DUA06060 

15 CONTINUE DUA06070 


************************************************************************00^06080 


** INITIAL CALCULATIONS DUA06090 

************** ********************************************************** 00^00 ^ 00 


DO 00100 I = 1 , PDTX, 1 
DTAARY(I) = DTXARX(I) 

00100 CONTINUE 

LAMBDA = C/FREQ 
K = 2*PI/LAMBDA 
MAXMX = PMXMX 
MAXMY = PMXMY 
MAXSX = PMXSX 
MAXSY = PMXSY 
MAGSR = 0 

DO 00200 V = 1,3,1 

R1TMP = REFXYZ (V) - FEDXYZ(V) 

SR(V) = R1TMP 

MAGSR = MAGSR + R1TMP*R1TMP 
00200 CONTINUE 

MAGSR=SQRT (MAGSR) 

DO 00300 V = 1,3,1 
USR(V) = SR(V) /MAGSR 
00300 CONTINUE 

******************* *************************************************** 


****** j^ipyy . FAR-FIELD LIMIT POINTS********************************** 
****** E-PLANE : 90 - 270 PHI CUTS 

****** H-PLANE : 0 - 180 PHI CUTS 

****** 45-PLANE: 45 - 225 PHI CUTS 

********************************************************************** 


BEG PH I = 74.99858*pi/180. 

ENDPHI = 254.99858*pi/180. 
IDXPHI = 1 
STPPHI = 1 

INCPHI = (ENDPHI-BEGPHI) /IDXPHI 


* INCPHI = 0 

********************************************************************** 


BEGTHE 

ENDTHE 

IDXTHE 

STPTHE 


0 . 

5.*PI/180 

50 

1 


INCTHE = (ENDTHE-BEGTHE) /IDXTHE 

*********************************************************************** 
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** Calculate Current Densities on the SubRef lector Resulting from SourceDUA06490 
************* *********************************************************** 00^00000 


INCSXY = INCSX*INCSY 

DUA06510 

XS = XSO - INCSX 

DUA06520 

DO 00700 I = 1, MAXSX, 1 

DUA06530 

XS = XS + INCSX 

DUA06540 

YS = YSO - INCSY 

DUA06550 

DO 00600 J = 1, MAXSY, 1 

DUA06560 

YS = YS + INCSY 

DUA06570 

ZS = SUBARY ( IDXZ, J, I) 

DUA06580 

MSKFAC = SUBARY (IDXMSK,J, I) 

DUA06590 

IF (MSKFAC .EQ. 1.) THEN 

DUA06600 
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00500 


00510 


00520 


00530 


00540 


00550 


MAGS I = 0. 

DO 00500 V = 1,3,1 

RITMP = SUBXYZ (V) - FEDXY2(V) 

SI (V) = RITMP 

MAGS I = MAGS I + R1TMP*R1TMP 
CONTINUE 

MAGSI=SQRT (MAGSI) 

IF (MAGSI .EQ. 0.) THEN . , 

STOP 

ENDIF 

DO 00510 V = 1,3,1 
USI(V) = SI (V) /MAGSI 
CONTINUE 

PVX = (-USIX*USIY) 

PVY = ( USIX*USIX + USIZ*USIZ) 

PVZ = (-USIY*USIZ) 

MAGPV = SQRT (PVX*PVX + PVY*PVY + PVZ*PVZ) 
DO 00520 V = 1,3,1 
UPV(V) = PV (V) /MAGPV 
CONTINUE 

COSTHE = D0T(USR,UPV,3) 

SCALE = FDPTRN(USR, USI , Q, COZ, FLG) /ETA 
IF ( (FLG .EQ. 1) 

.OR. (SCALE .EQ. 0.)) THEN 

STOP 


ENDIF 

CALL CROSS(HV,USI,UPV) 

MAGHV = SQRT (HVX*HVX + HVY*HVY + HVZ*HVZ) 
DO 00530 V = 1,3,1 
UHV(V) = HV(V) /MAGHV 
CONTINUE 

NX - SUBARY (IDXZX, J , I) 

NY = SUBARY ( IDXZY, J , I) 

NZ = -1 

MAGNRM = SQRT (NX*NX + NY*NY +1) 

DO 00540 V = 1,3,1 

RITMP = NORM(V) /MAGNRM 

NORM(V) = RITMP 

SCDARY (V , IDXNRM, J , I) = RITMP 
CONTINUE 

CALL CROSS (JV, NORM, UHV) 

MAGJV = SQRT (JVX*JVX + JVY*JVY + JVZ*JVZ) 
MAGJV =1. 

DO 00550 V = 1,3,1 

RITMP = JV (V) /MAGJV 

UJV(V) = RITMP 

SCDARY (V, IDXJVX, J , I) = RITMP 
CONTINUE 

IF (MAGSI .EQ. 0.) THEN 

STOP 

ENDIF 


SCDARY (IDXMNM, IDXDTX, J, 

SCDARY (IDXMSI, IDXDTX, J, 

SCDARY(IDXAOT, IDXDTX, J, 

SCALE 

COSKR 

SINKR 

JFLD(IDXVCX, IDXRMJ) = 
JFLD(IDXVCY, IDXRMJ) = 
JFLD(IDXVCZ, IDXRMJ) = 


I)= MAGNRM 
I)= MAGSI 
I)=(2.*SCALE) 
2.*SCALE/MAGSI 
COS (K*MAGSI) 

SIN(K*MAGSI) 

(SCALE*UJ V ( IDXVCX) *COSKR) 
(SCALE*UJV (IDXVCY)*COSKR) 
(SCALE*UJ V ( IDXVCZ) *COSKR) 
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JFLD(IDXVCX, IDXIMJ) = -(SCALE*UJV(IDXVCX)*SINKR) 

DUA07210 



JFLD(IDXVCY, IDXIMJ) = -(SCALE*UJV(IDXVCY)*SINKR) 

DUA07220 



JFLD(IDXVCZ, IDXIMJ) = - (SCALE*UJV(IDXVCZ)*SINKR) 

DUA07230 



ELSE 

DUA07240 


ENDIF 

DUA07250 

00599 

CONTINUE 

DUA07260 

00600 

CONTINUE 

DUA07270 

00700 CONTINUE 

DUA07280 

***i(******ic*i<**********'k**i<******'k*'k***i<*****ic*****itik****it*****i'********Q{)i\o'j290 

** Calculate 

Induced Magnetic Field on Main Reflector by Sub Reflector. 

DUA07300 

ic***********************************************************************^^^^ 

INCMXY 

= INCMX*INCMY 

DUA07320 

XM 


= XMO - INCMX 

DUA07330 

DO 

01900 IP = l.MAXMX.l 

DUA07340 


XM 

= XM + INCMX 

DUA07350 


YM 

= YMO - INCMY 

DUA07360 


DO 

01800 JP = l.MAXMY.l 

DUA07370 



YM = YM + INCMY 

DUA07380 



ZM = MAIARY (IDXZ.JP, IP) 

DUA07390 



MSKFAC = MAIARY (IDXMSK, JP, IP) 

DUA07400 



IF (MSKFAC .EQ. 1.) THEN 

DUA07410 



NX = - (MAIARY (IDXZX,JP, IP)) 

DUA07420 



NY = - (MAI ARY (IDXZY, JP, IP) ) 

DUA07430 



NZ = +1. 

DUA07440 



MAGNRM = SQRT (NX*NX + NY*NY +1.) 

DUA07450 



DO 01100 V = 1,3,1 

DUA07460 



NORM(V) = NORM (V) /MAGNRM 

DUA07470 



DO 01000 W = 1,2,1 

DUA07480 



HFLD(V, W) = 0. 

DUA07490 

01000 


CONTINUE 

DUA07500 

01100 


CONTINUE 

DUA07510 



XS = XSO - INCSX 

DUA07520 



DO 01700 I = 1 ,MAXSX, 1 

DUA07530 



XS = XS + INCSX 

DUA07540 



YS = YSO - INCSY 

DUA07550 



DO 01600 J = 1 ,MAXSY, 1 

DUA07560 



YS = YS + INCSY 

DUA07570 



ZS = SUBARY (IDXZ.J, I) 

DUA07580 



MSKFAC = SUBARY(IDXMSK, J , I) 

DUA07590 



IF (MSKFAC .EQ. 1.) THEN 

DUA07600 



DS = SCDARY(IDXMNM, IDXDTX, J, I)*INCSXY 

DUA07610 



MAGS I = SCDARY(IDXMSI, IDXDTX, J, I) 

DUA07620 



SCALE = SCDARY(IDXAOT, IDXDTX, J, I) 

DUA07630 



MAGR1 = 0. 

DUA07640 



DO 01200 V = 1,3,1 

DUA07650 



R1TMP = (MAIXYZ (V) - SUBXYZ(V)) 

DUA07660 



R1 ( V) = R1TMP 

DUA07670 



MAGR1 = MAGR1 + R1TMP*R1TMP 

DUA07680 

01200 


CONTINUE 

DUA07690 



MAGR1 = SQRT (MAGR1) 

DUA07700 



DO 01300 V = 1,3,1 

DUA07710 



UR1(V) = R1 (V) /MAGR1 

DUA07720 



UHV(V) = SCDARY (V, IDXJVX, J , I) 

DUA07730 

01300 


CONTINUE 

DUA07740 



KR =(K*(MAGR1 + MAGSI)) 

DUA07750 



RR = MAGR1*MAGSI 

DUA07760 



COSKR = +COS(KR)/RR 

DUA07770 



SINKR = -SIN(KR)/RR 

DUA07780 

* 


CALL CROSS (INTG, SCDARY ( IDXVCX , IDXJVX , J , I ) , UR1 ) 

DUA07790 



CALL CR0SS(INTG,UHV,UR1) 

DUA07800 
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MAG I NT = 1. 

MAGINT = SQRT ( ( (INTX*INTX) 
+(INTY*INTY) 
+(INTZ*INTZ))) 
INTX = INTX/MAGINT 


DUA07810 

DUA07820 

DUA07830 

DUA07840 

DUA07850 


INTY = INTY/MAGINT DUA07860 

INTZ = INTZ/MAGINT DUA07870 

DO 01500 V = 1,3,1 DUA07880 

DO 01400 W = 1,2,1 DUA07890 

0 HFLD(V,W) = HFLD(V,W) DUA07900 

1 +SCALE*PSI (W)*INTG(V)*DS DUA07910 

01400 CONTINUE DUA07920 

01500 CONTINUE DUA07930 

ENDIF DUA07940 

01599 CONTINUE DUA07950 

01600 CONTINUE DUA07960 

01700 CONTINUE DUA07970 

MCDARY (IDXVCX, IDXUNM, JP, IP) = NX DUA07980 

MCDARY(IDXVCY, IDXUNM.JP, IP) = NY DUA07990 

MCDARY (IDXVCZ, IDXUNM, JP, IP) = NZ DUA08000 

CALL SCALER (NORM, NORM, 2.) DUA08010 

CALL CROSS ( J FLD ( I DX VCX , IDXRMJ) , NORM,HFLD(IDXVCX, IDXRMJ) )DUA08020 
CALL CROSS (JFLD(IDXVCX, IDXIMJ) , NORM, HFLD (IDXVCX, IDXIMJ) )DUA08030 
DO 01720 V = 1,3,1 DUA08040 

DO 01710 W = 1,2,1 DUA08050 

MCDARY (V, W, JP, IP) = JFLD(V,W) DUA08060 

01710 CONTINUE DUA08070 

01720 CONTINUE DUA08080 

******************* ****************************************** ***********Q(j^QgQgQ 

ELSE DUA08100 

JFLD(IDXVCX, IDXRMJ) = 0. DUA08110 

JFLD(IDXVCY, IDXRMJ) = 0. DUA08120 

JFLD(IDXVCZ, IDXRMJ) = 0. DUA08130 

JFLD(IDXVCX, IDXIMJ) = 0. DUA08140 

JFLD(IDXVCY, IDXIMJ) = 0. DUA08150 

JFLD(IDXVCZ, IDXIMJ) = 0. DUA08160 

MAGJV = 0. DUA08170 

ENDIF DUA08180 

01799 CONTINUE SuA08200 

01800 CONTINUE DUA08210 

01900 CONTINUE DUA08220 

****************************************** ******************************py^gg2 ^ Q 

*** FAR FIELD ANTENNA PATTERN COMPUTATION DUA08240 

*********************************************************************** iAnQ9t;n 

PTTMIN = +1.E+38 DUA08260 

PTTMAX = -l.E+38 DUA08270 

ANGPHI = BEGPHI - INCPHI DUA08280 

DO 02600 ANGPHX = 0, IDXPHI .STPPHI DUA08290 

ANGPHI = ANGPHI + INCPHI DUA08300 

SINPHI = SIN (ANGPHI ) DUA08310 

COSPHI = COS(ANGPHI) DUA08320 

ANGTHE = BEGTHE - INCTHE DUA08330 

DO 02500 ANGTHX = 0, IDXTHE.STPTHE DUA08340 

ANGTHE = ANGTHE + INCTHE DUA08350 

SINTHE = SIN(ANGTHE) DUA08360 

COSTHE = COS (ANGTHE) DUA08370 

RFFX = SINTHE*COSPHI DUA08380 

RFFY = SINTHE*SINPHI DUA08390 

RFFZ = COSTHE DUA08400 
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oo ro i— ‘ wi\)m wno 


SUM (1,1) = 0. 

SUM(1 , 2) = 0. 

SUM(2, 1) = 0. 

SUM(2,2) = 0. 

SUM(3, 1) = 0. 

SUM(3,2) = 0. 

XM = XM0 - INCMX 
DO 02400 IP = 1 , MAXMX, 1 
XM = XM + INCMX 
YM = YM0 - INCMY 
DO 02300 JP = 1 ,MAXMY, 1 
YM = YM + INCMY 
ZM = MAI ARY (IDXZ, JP, IP) 

KR = (K*( (RFFX*XM)+(RFFY*YM)+(RFFZ*ZM)) ) 

COSKR = COS(KR) 

SINKR = SIN (KR) 

MSKFAC = MAIARY (IDXMSK.JP, IP) 

IF(MSKFAC.EQ.O)GO TO 2300 

NORM (IDXVCX) = -MAIARY (IDXZX.JP, IP) 

NORM (IDXVCY) = -MAIARY(IDXZY, JP, IP) 
NORM(IDXVCZ) = +1. 

DS = INCMXY*SQRT ( NORM ( I DXVCX ) *NORM( IDXVCX) 
+NORM(IDXVCY)*NORM(IDXVCY) 
+NORM(IDXVCZ) *NORM(IDXVCZ) ) 

CMPTMP (IDXRMJ) = 

RFFX*MCDARY (IDXVCX, IDXRMJ, JP, IP) 
+ RFFY*MCDARY (IDXVCY , IDXRMJ , JP , I P) 
+ RFFZ*MCDARY(IDXVCZ, IDXRMJ, JP, IP) 
CMPTMP(IDXIMJ) = 

RFFX*MCDARY ( IDXVCX , IDXIMJ , JP , I P) 
+ RFFY*MCDARY (IDXVCY, IDXIMJ, JP, IP) 
+ RFFZ*MCDARY(IDXVCZ, IDXIMJ.JP, IP) 
TMRX = CMPTMP ( IDXRMJ )*RFFX 
TMRY = CMPTMP ( IDXRMJ )*RFFY 
TMRZ = CMPTMP ( IDXRMJ ) *RFFZ 
TMIX = CMPTMP ( IDXIMJ) *RFFX 
TMIY = CMPTMP( IDXIMJ) *RFFY 
TMIZ = CMPTMP( IDXIMJ) *RFFZ 
TMRX = MCDARY (IDXVCX, IDXRMJ, JP, IP) - TMRX 
TMRY = MCDARY(IDXVCY, IDXRMJ, JP, IP) - TMRY 
TMRZ = MCDARY (IDXVCZ, IDXRMJ, JP, IP) - TMRZ 
TMIX = MCDARY(IDXVCX, IDXIMJ, JP, IP) - TMIX 
TMIY = MCDARY(IDXVCY, IDXIMJ.JP, IP) - TMIY 
TMIZ = MCDARY(IDXVCZ, IDXIMJ.JP, IP) - TMIZ 
TTMRX = ( (TMRX*COSKR) - (TMIX*SINKR) ) 

TTMRY = ( (TMRY*COSKR) - (TMIY*SINKR) ) 

TTMRZ = ( (TMRZ*COSKR) - (TMIZ*SINKR) ) 

TTMIX = ( (TMRX*SINKR)+(TMIX*COSKR)) 

TTMIY = ( (TMRY*SINKR)+(TMIY*COSKR)) 

TTMIZ = ( (TMRZ*SINKR)+(TMIZ*COSKR) ) 

JVRX = TTMRX*DS 
JVRY = TTMRY*DS 
JVRZ = TTMRZ*DS 
JVIX = TTMIX*DS 
JVIY = TTMIY*DS 
JVIZ = TTMIZ*DS 
SUMRX = SUMRX + JVRX 
SUMRY = SUMRY + JVRY 
SUMRZ = SUMRZ + JVRZ 
SUMIX = SUMIX + JVIX 
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02299 

02300 
02400 

0 

1 

2 

0 

1 

2 

0 

1 

0 

1 


SUMIY = SUMIY + JVIY 
SUMIZ = SUMIZ + JVIZ 

CONTINUE 

CONTINUE 

CONTINUE 

ETHRE = SUMRX*COSTHE*COSPHI 
+ SUMRY*COSTHE*SINPHI 

- SUMRZ*SINTHE 

ETHIM = SUMIX*COSTHE*COSPHI 
+ SUMIY*COSTHE*SINPHI 

- SUMIZ*SINTHE 
EPHRE = - SUMRX*SINPHI 

+ SUMRY*COSPHI 
EPHIM = - SUMIX*SINPHI 
+ SUMIY*COSPHI 


*********************************************************************** 


**********£j|^^j^ AND £pnj COMPUTED ************************************ 
********** T rAnsf 0 RM TO CO AND CO-POL LUDWIG’S DEFINITION*************** 
*********************************************************************** 

AREFR=SINPHI*ETHRE+COSPHI*EPHRE 

AREFI=SINPHI*ETHIM+COSPHI*EPHIM 

ACRPR=-COSPHI*ETHRE+SINPHI*EPHRE 

ACRPI=-COSPHI*ETHIM+SINPHI*EPHIM 

********************************************************************** 

AMGREF=AREFR**2+AREFI**2 

AMGCRP=ACRPR**2+ACRPI**2 

********** C 0 AND CROSS POL FIELDS ARE COMPLETED*********************** 
********************************************************************** 

*************************** ****jotal field**************************** 

* 0 R1TMP = ETHRE*ETHRE + ETHIM*ETHIM 

* 1 + EPHRE*EPHRE + EPHIM*EPHIM 

***************PLqj THE CO POL FIELDS********************************* 

R1TMP = AMGREF 

***************PI_qj y^|E CROSS POL FIELDS****************************** 

* R1TMP = AMGCRP 

********************************************************************** 
PTTRN(IDXPWR,ANGTHX, ANGPHX) = R1TMP 
IF (R1TMP .LT. PTTMIN) THEN 
PTTMIN = R1TMP 
ENDIF 

IF (R1TMP .GT. PTTMAX) THEN 
PTTMAX = R1TMP 
ENDIF 

02500 CONTINUE 
02600 CONTINUE 

PRAD = ( (2. *PI) / (ETA* ( (2. *Q)+1) ) ) 

RIFCTR = ( ( (K*K*K*K) *ETA) / (256. * (PI*PI*PI*PI) ) ) 

RINTNS = (RIFCTR*PTTMAX) 

DIRCTV = ( (4. *PI*RINTNS) / (PRAD) ) 

DIRCTV = (10*AL0G10 (DIRCTV) ) 

WRITE (19 ,897) DIRCTV 
897 FORMAT (5X, F15.5) 

ADBMNX = 10. *AL0G10( PTTMIN) 

ADBMXX = 10 . *AL0G10 (PTTMAX) 

RDBMNX - (ADBMNX - ADBMXX) 

ANGPHI = BEGPHI - INCPHI 
DO 02800 ANGPHX = 0, IDXPHI , STPPHI 
ANGPHI = ANGPHI + INCPHI 
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ANGTHE = BEGTHE - INCTHE 
DO 02700 ANGTHX = 0, IDXTHE, STPTHE 
ANGTHE = ANGTHE + INCTHE 

R1TMP = 10.*ALOG10( (PTTRN (I DXPWR, ANGTHX .ANGPHX) ) ) 
PTTRN(IDXADB, ANGTHX, ANGPHX) = (R1TMP ) 

PTTRN ( IDXRDB , ANGTHX , ANGPHX) = (R1TMP - ADBMXX) 


=DUA09670 


*******Q£^£P^-j-^ OUTPUT ARRAYS***************************************DUA09680 
************************************************************************DU^09690 


895 FORMAT ( 

02700 CONTINU 

02800 CONTINUE 
DTXARY ( 2) 
DTXARY ( 3) 
DTXARY ( 4) 
DTXARY ( 5) 
DTXARY ( 6) 
DTXARY ( 7) 
DTXARY ( 8) 
DTXARY ( 9) 
DTXARY (10) 
DTXARY (11) 
DTXARY (12) 


ZETA=ANGTHE*180/PI 

APHI=ANGPHI*180/PI 

WRITE ( 19 , 895) PTTRN ( IDXRDB, ANGTHX , ANGPHX) , ZETA, APHI 
FORMAT (5X,F15.8,3X,2(F15.8,3X)) 

CONTINUE 


ENDPHI 

REAL(IDXPHI) 

REAL(STPPHI) 

INCPHI 

BEGTHE 

ENDTHE 

REAL(IDXTHE) 

REAL(STPTHE) 

INCTHE 

PTTMIN 

PTTMAX 
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END DUA09870 

************************************************************ * ********** QU AO 9880 
**********r E A|_ FUNCTION FDPTRN RETURN FEED PATTERN********************* DUA09890 
**************************** ********************************** ********* DUA09900 
REAL FUNCTION FDPTRN (THETA, PHI, RHO,COZ,ERR) DUA09910 

REAL THETA (3) DUA09920 

REAL PHI (3) DUA09930 

REAL RHO DUA09940 

REAL COZ DUA09950 

INTEGER ERR DUA09960 

REAL DOT DUA09970 

EXTERNAL DOT DUA09980 

REAL DOTVAL DUA09990 

DOTVAL = D0T(THETA,PHI,3) DUA10000 

IF (DOTVAL .LT. 0.) THEN DUA10010 

ERR = 1 DUA10020 

ERR = 1 DUA10030 

COZ =0. DUA10040 

FDPTRN = 0. DUA10050 

ELSE DUA10060 

ERR = 0 DUA10070 

COZ = DOTVAL DUA10080 

FDPTRN = (DOTVAL) **RHO DUA10090 

ENDIF DUA10100 

RETURN DUA101 10 

END DUA10120 

************************************************************************DUAioi30 
** REAL FUNCTION DOT() ! Returns Real Value of DOT PRODUCT A and B DUA10140 
******************************************************************** ****qu^ 1 50 
REAL FUNCTION DOT (A,B,N) DUA10160 

INTEGER N DUA10170 

REAL A(N) DUA10180 

REAL B(N) DUA10190 

INTEGER I DUA10200 


ERR 

COZ 

FDPTRN 

ELSE 

ERR 

COZ 

FDPTRN 

ENDIF 

RETURN 


******************************************************************** 

** REAL FUNCTION DOT() ! Returns Real Value of DOT PRODUCT A and 
******************************************************************** 


0 

DOTVAL 

(DOTVAL) **RHO 


REAL FUNCTION DOT (A,B,N) 
INTEGER N 
REAL A(N) 

REAL B(N) 

INTEGER I 
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REAL SUM DUA10210 

SUM = 0. DUA10220 

DO 00100 I = 1 , N, 1 DUA10230 

SUM = SUM + A(I)*B(I) DUA10240 

00100 CONTINUE DUA10250 

DOT = SUM DUA10260 

RETURN DUA10270 

END DUA10280 


************************************************************************01 JA10290 


** SUBROUTINE CROSS ! Performs C = AxB DUA10300 


************************************************************************00/^03 1 0 


SUBROUTINE CROSS(C,A,B) 

DUA10320 

REAL 

C (3) 

DUA10330 

REAL 

A ( 3 ) 

DUA10340 

REAL 

B(3) 

DUA10350 

c(i) 

= +((A(2)*B(3))-(A(3)*B(2))) 

DUA10360 

C(2) 

- -((A(1)*B(3))-(A(3)*B(1))) 

DUA10370 

C ( 3 ) 

= +((A(1)*B(2))-(A(2)*B(1))) 

DUA10380 

END 


DUA10390 

********************** **************************************************00^2 0400 

** SUBROUTINE 

SCALER ! Performs C = A*<SCALER> 

DUA10410 

****** ************************************************************ ******00/\ 10420 

SUBROUTINE SCALER(C, A, SCALEX) 

DUA10430 

REAL 

C(3) 

DUA10440 

REAL 

A (3) 

DUA10450 

REAL 

SCALEX 

DUA10460 

c(i) 

= SCALEX*(A(1) ) 

DUA10470 

C (2) 

= SCALEX* (A(2) ) 

DUA10480 

C ( 3 ) 

= SCALEX* (A(3) ) 

DUA10490 

RETURN 


DUA10500 

END 


DUA10510 

********************************************************************** **00^ ]^Q520 

** SUBROUTINE 

VECADD ! Performs C = A+B 

DUA10530 


***************************************************************** ******* 00 ^ | 0540 


SUBROUTINE VECADD(C, A, B) 

DUA10550 

REAL 

C (3) 

DUA10560 

REAL 

A ( 3 ) 

DUA10570 

REAL 

B ( 3 ) 

DUA10580 

C(l) 

= (A(l)-B(l)) 

DUA10590 

C(2) 

= (A(2)+B(2)) 

DUA10600 

C (3) 

- (A(3)+B(3) ) 

DUA10610 

RETURN 


DUA10620 

END 


DUA10630 

************************************************************************ 00^20540 

** SUBROUTINE 

VECSUB ! Performs C = A-B 

DUA10650 


*********************************************************************** *Q 0 /\ \ 0660 


SUBROUTINE VECSUB(C. A, B) 

DUA10670 

REAL 

C (3) 

DUA10680 

REAL 

A (3) 

DUA10690 

REAL 

B (3 ) 

DUA10700 

C(l) 

= (A(l)-B(l)) 

DUA 10710 

C(2) 

- (A (2) -B(2) ) 

DUA10720 

C (3) 

= (A (3) -B (3) ) 

DUA10730 

RETURN 


DUA10740 

END 


DUA10750 
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PROGRAM FFPLOT 

DIMENSION X (10000) , Y (10000) , VARS (20) 

DIMENSION XPL(IOOO) , YPL ( 1000) 

CHARACTERS CH/ 1 DIRECTIVITY =7 
CHARACTER*2 ADB/'DBV 
CHARACTER*5 DIR(l) 

C*****THIS PROGRAM CAN BE USED TO PLOT THE ANTENNA FAR-FIELD PATTERN 
C***** (E-PLANE OR H-PLANE CUTS) 
c *****!axis , NUM , Y , RTNARR : PARAMETERS IN SCLBK2 
C ***** IVAR; PARAMETERS IN GPL0T3 
INTEGER * 4 IAXIS /O/ 

INTEGER * 4 IVARS (20) 

INTEGER * 4 NUM/10000/ 

CHARACTERS XTITLE(5) / ' ELEV ' , ' ATIO ' , ' N AN'.'GLE ' , 'DEG. ' / 
CHARACTERS YTITLE (6)/' RELA', 'TIVE', ' AMP' , 1 LITU' , 'DE ' , ' (DB) 7 
C*****NP : TOTAL NO. OF POINTS ; ZM: MAXM. VIEWING ANGLE(DEG. ) 

READ(19, 756)DDIR 
756 FORMAT (5X, F15.5) 

****** 2 (NFF+1) 

NP=102 

ZM=2. 


DO 15 J=1 , NP 

C*****X : ANGLE POSITIONS(DEG. ) ; Y: RELATIVE FAR FLD. AMPLITUDES (DB) 
READ(19, 300) YPL(J) , XPL ( J ) ,DUM1 
300 FORMAT (5X,F15.8,3X,2(F15.8,3X)) 

15 CONTINUE 

DO 98 J=1 , 51 
Y (J) = (YPL(103-J) ) 

X(J) =— XPL ( 1 03 - J ) 

Y (J+51) = (YPL(J) ) 

X(J+51) =XPL(J) 

98 CONTINUE 

C ***** SCL K K 2 : GRAPH3D ROUTINE TO FIND MIN, MAX IN DATA 
C***** o . y-COORDINATE ; NUM : DIMENSION OF Y-ARRAY ; Y : Y-ARRAY 
C*****RTNARR(2) : DIMENSION TO STORE Y(MIN) , Y(MAX) VALUES 
£*****PE^pp^g[: j|_|£ FIELD VALUES*********************************** 
C*****U XTRM :GRAPH3D ROUTINE; DEFINES EXTREME POSITIONS OF A 3D PLOT 
C*****8 ; TOTAL NO. OF VARIABLES ; 0 : CARTESIAN ; (-ZM.ZM) :(XMIN f XMAX) 
C*****(-80.,0) : (YMIN, YMAX) ; (0..0.) :(ZMIN,ZMAX) 

CALL UXTRM(8, 0,-ZM,ZM, -54. ,0., 0.0,0.) 

C*****UMAPF :GRAPH3D ROUTINE , DEFINE MAPPING TO TRANSFORM FROM USER 
C*****TO RELATIVE UNITS. 

C***** o CARTESIAN ; 1. :ONE VARS, DEFAULT ; 0 :N0 LOG SCALE 
CALL UMAPF (0, 1 . , 0) 

C*****XA XIS 3 : GRAPH3D ROUTINE , DEFINES X-AXIX COORDINATE 
C***** VARS: 1=T0TAL NO OF VARS ; 2 3 4=X1 Y1 Z1 ; 5 6 7=X2 Y2 12 
C***** 8=USER UNIT (1 .) ; 9=N0. OF INTERVALS ; 10=GRID 0PTI0N(1.) 

C***** 11=DRAW PARALLEL TO Y-AXIS ; 12=VARS(9)+1 ; 13=SIZE OF LABEL 
C***** 14= (DIR. OF X AXIS) (CENTERED AT GRID) (CLOCKWISE TO AXIS) 

C***** 15=AXIS SETTING IS NOT COMPLETE 
VARS(l) =15 
VARS(2) =-ZM 
VARS(3)=-54. 

VARS(4) =0.0 
VARS(5)=ZM 
VARS (6) =-54. 

VARS(7)=0. 


FFP00010 

FFP00020 

FFP00030 

FFP00040 

FFP00050 

FFP00060 

FFP00070 

FFP00080 

FFP00090 

FFP00100 

FFP00110 

FFP00120 

FFP00130 

FFP00140 

FFP00150 

FFP00160 

FFP00170 

FFP00180 

FFP00190 

FFP00200 

FFP00210 

FFP00220 

FFP00230 

FFP00240 

FFP00250 

FFP00260 

FFP00270 

FFP00280 

FFP00290 

FFP00300 

FFP00310 

FFP00320 

FFP00330 

FFP00340 

FFP00350 

FFP00360 

FFP00370 

FFP00380 

FFP00390 

FFP00400 

FFP00410 

FFP00420 

FFP00430 

FFP00440 

FFP00450 

FFP00460 

FFP00470 

FFP00480 

FFP00490 

FFP00500 

FFP00510 

FFP00520 

FFP00530 

FFP00540 

FFP00550 

FFP00560 

FFP00570 

FFP00580 

FFP00590 
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VARS (8) =1. FFP00600 

VARS(9)=10. FFP00610 

VARS(IO) =1 . FFP00620 

VARS (11) =4. FFP00630 

VARS(12) =11 . FFP00640 

VARS(13)=20. FFP00650 

VARS (14) =211 . FFP00660 

VARS(15)=0. FFP00670 

CALL XAXIS3 (VARS) FFP00680 

C ***** VARS ; 2 3 4 =X1 Y1 Z 1 ; 5 6 7=X 2 Y2 Z2 ; 8=USER UNIT ; 9=NO. OF FFP00690 
C ***** INTERVALS . io=GRID OPTION ; 11=DRAW PARALLEL TO X-AXIS ; 12= FFP00700 

C ***** VARS ( 9 ) +1 . 13=SIZE OF LABEL ; 14=DIR. OF X-AXIS ; 15=AXIS SETTINGFFP00710 
C*****COMPLETE. FFP00720 

VARS(2) =-ZM FFP00730 

VARS(3)=-54. FFP00740 

VARS(4) =0. FFP00750 

VARS(5) =-ZM FFP00760 

VARS(6) =0. FFP00770 

VARS(7)=0. FFP00780 

VARS(8)=1. FFP00790 

VARS(9) =9. FFP00800 

VARS(10)=1. FFP00810 

VARS ( 1 1 ) =3 . FFP00820 

VARS (12)=10. FFP00830 

VARS(13) =20. FFP00840 

VARS (14)=212. FFP00850 

VARS(15) =1 . FFP00860 

CALL YAXIS3 (VARS) FFP00870 

C ***** TITL E3 : GRAPH3D ROUTINE ; PRINTS TITLE OF X-AXIS FFP00880 

C***** 4=X-AXIS ; 20=X- ALPHANUMERIC DIMENSION ; 15=CHARACTER SIZE FFP00890 

CALL TITLE3 (4, 20, 15, XTITLE, 0. , 1 . , 0. ) FFP00900 

C ***** TITLE3 : GRAPH3D ROUTINE ; PRINTS Y-AXIS TITLE FFP00910 

C***** 3=Y-AXIS ; 18=Y-ALPHANUMERIC DIMENSION ; 15=CHARACTER SIZE FFP00920 

CALL TITLE3 (3,24,15,YTITLE,-1.,0.,0.) FFP00930 

C ***** GPL0T3 . GRAP H3D ROUTINE ; TO PLOT A CURVE WITH POINT OR VECTOR FFP00940 

C ***** IVARS . DIMENSION OF IVARS ; 2=N0. OF POINTS, EXACT ; 3=N0 Z-AXISFFP00950 
4=D0 N oT CALL AXIS ROUTINES ; 5=P0INT PLOT ; 6=SYMB0L FREQUENCY FFP00960 


C***** 7=SIZE OF SYMBOL ; 8=EXACT MIN-MAX INTERVAL 
IVARS (1) =8 
IVARS(2) =NP 
IVARS(3) =0 
IVARS(4) =0 
IVARS(5) =0 
IVARS (6) =1 
IVARS(7)=15 
IVARS (8) =1 

CALL GPL0T3(IVARS,X, Y) 

CALL CHARS3 (13 , CH, 3. ,10. 5,0. ,25,1.) 

CALL NUMBER (4, DDIR, 5, 2, DIR) 

CALL CHARS3(5,DIR,7. ,10. 5,0., 25,1.) 

CALL CHARS3(2,ADB,9. ,10.5,0. ,25,1.) 

C ***** GVI EW : GRAPH3D ROUTINE ; IDENTIFIES VIEWING ENVIRONMENT 
C***** 1=DEFAULT VALUES FOR THREE REMAINING VARIABLES 
CALL GVIEW(l) 

C ***** WINDW . GRAPH3D ROUTINE ; SPECIFY DIMENSION OF VIEW WINDOW 
C***** 6=T0TAL NO OF VARIABLES ; 0=L0WEST OF THE PARAMETER RANGE 

C ***** UMIN=MIN> VALUE 0F N0< 0F RELATIVE UNITS FR0M VIEW REFEREN ce 

C*****UM AX =MAX. VALUE OF NO. OF RELATIVE UNITS FROM VIEW REFERENCE 

Q*****y^ j » ii ii m ii H ii n it ii 

Q*****^|^^^ = 'I II II II II M II II II " 


FFP01180 

FFP01190 
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CALL WINDW(6 t 0, -6.5,6.5,-6.5,6.5) FFP01200 

C ***** DISPLA : GRAPH3D ROUTINE ; DISPLAYS INTERNAL BUFFER FFP01210 

C***** l=OPTION TO CLEAR BUFFER FFP01220 

CALL DISPLA ( 1 ) FFP01230 

C*****TERM : GRAPH3D ROUTINE ; REQUIRED TO CLOSE THE GRAPHICS FFP01240 

CALL TERM FFP01250 

STOP FFP01260 

END FFP01270 
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/* EXEC DUAL REFLECTOR */ 
"GRAPH3D" 

SETUP FTN 

"FI 19 DISK DUALREF 0UT19 Al" 
"FI 15 DISK DUALREF OUT15 Al" 
"FI 16 DISK DUALREF OUT16 Al" 
"FI 17 DISK DUALREF OUT17 Al" 
"FI 18 DISK DUALREF 0UT18 Al" 
"LOAD DRSG (CLEAR START" 

"LOAD DUALREF (CLEAR START" 
"LOAD FFPLOT (CLEAR START" 
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DUAL REFLECTOR 
GEOMETRY 



Figure 1, Dual reflector configuration 
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Main Reflector 


XM 



Figure 2 , Generalized dual reflector geometry 
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Figure 3, Eulerian angles 
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RELATIVE AMPLITUDE (OB) 



Figure 4a, H-plane far-field antenna pattern 
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RELATIVE AMPIITUOE <OB) 



Figure 4b, E-plane far-field antenna pattern 
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